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I. INTRODUCTION 



What can one learn from measurements with the Sudbury Neutrino Observatory 
(SNO) What are the most likely quantitative results for each of the different experi- 
ments that can be carried out with SNO? The main goal of this paper is to help answer these 
questions by providing quantitative predictions for the most important diagnostic tests of 
neutrino oscillations that can be performed with SNO. 

SNO is not an experiment. Like LEP and Super-Kamiokande, SNO is a series of exper- 
iments. We calculate the currently-favored range of predictions for 10 quantities that are 
affected by neutrino oscillations and which SNO will measure. For the impatient reader, we 
list here the quantities that are sensitive to neutrino oscillations which we investigate (defi- 
nitions are given later in the text): first and second moments of the recoil energy spectrum, 
the charged current (CC), the neutral current (NC), and the neutrino-electron scattering 
rates, the difference between the day and the night rates for both the CC and the NC, the 
difference in the winter-summer CC rates, the neutral current (NC) to charged current (CC) 
double ratio, and the neutrino-electron scattering to CC double ratio. 

The simultaneous analysis of all the SNO results, measured values and upper limits, will 
be a powerful technique for constraining neutrino oscillation parameters. As an initial step 
in this direction, we analyze the combined results for five especially informative pairs of 
oscillation parameters. 

A. SNO reactions 

The SNO collaboration will study charged current (CC) neutrino absorption by deu- 
terium, 

v e + d — > p + p + e~ , (1) 
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neutral current (NC) neutrino disassociation of deuterium, 

v x + d — ► n +p + u' x , (x = e,/i, r), (2) 

and neutrino-electron scattering (ES), 

v x + e~ — > v' x + e~ , (x = e, /i, r). (3) 

The energy of the recoil electrons can be measured for the CC reaction, Eq. ([!]), and also for 
the ES reaction, Eq. For both these reactions, the operating energy threshold for the 
recoil electrons may be of order 5 MeV. The threshold for the NC reaction, Eq. @, is 2.225 
MeV. Just as for radiochemical solar neutrino experiments, there is no energy discrimination 
for the NC reaction. 

The Kamiokande || and Super-Kamiokande experiments || have performed precision 
studies of solar neutrinos using the neutrino-electron scattering reaction, Eq. (|3]). SNO will 
be the first detector to measure electron recoil energies as a result of neutrino-absorption, 
Eq. (fl|). We have presented in Ref. || detailed predictions of what may be observed with 
SNO for the CC (absorption) reaction. 

If there are no neutrino oscillations, i.e., 4>(v e ) = 0(total), then the ratios of the event 
rates in the SNO detector are calculated to be approximately in the following proportions: 
CC:NC:ES = 2.05:1.00:0.19, i.e., the number of CC events is expected to exceed the number 
of v-e scattering events by about a factor of 11. Since the NC efficiency is likely to be 
only about a half of either the CC or the ES efficiency [|IJ and currently favored oscillation 
solutions give <f>(v e ) ~ 0(total), the observed ratio of events in the SNO detector may actually 
be reasonably close to: CC:NC:ES ~ 2.0 : 0.5 : 0.2. 

In thinking about what SNO can do, it is useful to have in mind some estimated event 
rates for a year of operation. The Super-Kamiokande event rate [[| for neutrino electron 
scattering is 0.475 times the event rate that is predicted by the standard solar model ||. 
If there are no neutrino oscillations and the total solar neutrino flux arrives at earth in the 
form of u e with a 8 B neutrino flux of 0.475 times the standard model flux, then one expects 
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about 4.4 x 10 3 CC events per year in SNO above a 5 MeV threshold and about 1.1 x 10 3 NC 
events, while there should only be about 415 ES events. The above rates were calculated for 
a 5 MeV CC and ES energy threshold and for a 50% detection efficiency for NC events. For 
an 8 MeV threshold, the estimated CC rate is about 45% of the rate for a 5 MeV threshold 
and the ES rate is only about 28% of the 5 MeV threshold rate. For the currently favored 
oscillation solutions, the expected CC rates are typically of order 80% of the rates cited 
above and the NC rates are about a factor of two or three higher. 

B. What do we calculate? 

In this paper, we calculate the likely range of quantities that are measurable with SNO 
using a representative sample of neutrino parameters from each of the six currently allowed 
99% C.L. domains of two-flavor neutrino oscillation solutions. In other words, we explore 
what can be learned with SNO, assuming the correctness of one of the six neutrino oscil- 
lation solutions that is globally consistent at 99% C.L. with all of the solar neutrino 
experiments performed so far (chlorine J7], Kamiokande [fj], Super-Kamiokande 0, Sage ||, 
and GALLEX §). 

Table | lists the mixing angles and differences of mass squared for the six global best-fit 
solutions. Figure [l] and Figure | show the survival probabilities of the best-fit solutions as 

TABLE I. Best-fit global oscillation parame- 
ters. The differences of the squared masses are given 
in eV 2 . The survival probabilities that correspond to 
these best-fit solutions are shown in Fig. |l] and Fig. ^. 
Results are taken from Ref. Q 



Scenario 


Am 2 




sin 2 (26) 


LMA 


2.7 x 10" 


-5 


7.9 x lO" 1 


SMA 


5.0 x 10~ 


-6 


7.2 x 10~ 3 


LOW 


1.0 x 10~ 


-7 


9.1 x 10" 1 
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6.5 x 10" 


-11 


7.2 x 10" 1 
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4.4 x 10" 
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6.6 x 10~ 3 
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FIG. 1. The survival probabilities as a function of neutrino energy. For each of the six best-fit 
globally acceptable neutrino oscillation solutions discussed in Ref . Q , the figure shows the survival 
probability for electron type neutrinos as a function of energy. The results are averaged over one 
year. 



For each measurable quantity i, we express our predictions based upon neutrino oscilla- 
tion models in terms of the value predicted by an oscillation scenario divided by the value 
predicted by the combined standard electroweak model and the standard solar model. Thus 
for each measured quantity, i (like CC or NC event rate), we evaluate the expected range of 



*For the MSW solutions, there are small but perceptible differences in the computed survival 
probabilities which depend upon the neutrino production probability as a function of solar radius. 
The survival probabilities shown in Fig. [I] were computed by averaging the survival probability 
over the 8 B production region in the BP98 model ||. In order to portray more accurately the 
behavior at low energy, the survival probabilities for Fig. Q were computed by averaging over the 
p — p production region in the BP98 model. 
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the reduced quantity [i] 

_ (Observed Value); 
^ = (Standard Model Value);' ^ ' 

The reduced quantity [i] is by construction independent of the absolute value of the solar 
neutrino flux, which is used in calculating both the numerator and denominator of Eq. (|]) . 



t 1 1 1 1 1 — i — i — | 1 1 1 1 1 1 — i — i — |" 
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FIG. 2. Comparison of survival probabilities. For four of the best-fit oscillation solutions, the 
survival probabilities are compared using a linear energy scale. The differences in the energy 
dependence of the survival probability between the high energy region, ~ 8 MeV, and the low 
energy region, < 1 MeV, can be seen clearly on this figure. The MSW Sterile solution has an 
energy dependence similar to the SMA solution. For clarity, we have omitted the VACl solution. 
The parameters of the best-fit solutions are given in Table |. 

What fluxes are used in calculating the predicted rates (e.g., for charged current or 
electron-neutrino interactions) implied by different neutrino scenarios? We determine the 
best-fit ratio of the observed neutrino flux to the standard model flux by fitting to the 
Super-Kamiokande rate and observed recoil electron spectrum. The procedure is described 
in Sec. |T| following Eq. ([12]) and Eq. fll3|) (see especially the definition of f( 8 B)). We 
have not included explicit uncertainties in determining f{ 8 B) for a given pair of neutrino 
variables, Am 2 and sin 2 29, but we instead have allowed f( s B) to range over the wide set 
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of values obtained by applying our best-fit procedure at each point in the currently-allowed 
neutrino-parameter space. 

As we shall see, the most powerful diagnostics of neutrino oscillations are formed by 
considering the reduced double ratio of two measurable quantities, i and j, as follows: 

[i] _ (Observed Value)j/(Standard Model Value); 
[j] (Observed Value)j/ (Standard Model Value)j 

For example, the reduced double ratio of NC to CC rates is not only independent of the ab- 
solute flux of the solar neutrinos but is also insensitive to some experimental and theoretical 
uncertainties that are important in interpreting the separate [NC] and [CC] rates. 

We describe how we evaluate the uncertainties in Sec. |J. All of the calculated departures 
from the standard model expectations are small except for the double ratio of NC to CC, 
[NC]/[CC]. Therefore, the theoretical and the experimental uncertainties are important. 

We present in Sec. [Ill] the results predicted by the six oscillation solutions for the first and 
second moments of the shape of the CC recoil electron energy distribution. We summarize 
in Sec. [TV] the principal predictions for the CC rate, in Sec. |V| the predictions for the neutral 
current rate, and in Sec. [VT| the predictions for the neutrino-electron scattering rate. We 
then calculate the detailed predictions of the most important double ratios, the NC to 



charged current ratio, [NC]/[CC], in Sec. |VII| and the neutrino-electron scattering to CC 
ratio, [ES]/[CC], in Sec. [VII]. 



Up to this point in the paper, i.e., through Sec. [VII]] , we only discuss time-averaged 
quantities. In Sec. |X|, we present the predictions for the CC of the difference between the 
event rate observed at night and the event rate observed during the day. For the NC rate, 
there is also a small difference predicted between the night rate and the day rate if the MSW 
Sterile solution is correct. We analyze in Sec. [X]the seasonal effects in the CC rate. 

Section [XI] is a pairwise exploration of the discriminatory power gained by analyzing 
simultaneously the predictions and the observations of different smoking-gun indicators of 
neutrino oscillations. We consider in this section the joint analysis of variables like [NC]/[CC] 
versus the first moment of the CC energy spectrum, the day-night difference, or the neutrino- 



electron scattering rate. In Sec. |XII] , we summarize and discuss our principal conclusions. 
Since we evaluate so many different effects, we give in Sec. [Xill| our personal list of our top 
four conclusions. 



C. How should this paper be read? 

We recommend that the reader begin by looking at the figures, which give a feeling for 
the variety and the size of the various quantities that can be measured with SNO. Then we 
suggest that the reader jump directly to the end of the paper. The main results of the paper 



are presented in this concluding section; the summary given in Sec. |XI1] can be used as a 
menu to guide the reader to the detailed analyses that are of greatest interest to him or to 
her. 

This is the fifth in a series of papers that we have written on the potential of the Sudbury 
Neutrino Observatory for determining the properties of neutrino oscillations. The reader 
interested in details of the analysis may wish to consult these earlier works P,[T0|-|T!^ , which 
also provide a historical perspective from which the robustness of the predictions can be 
judged. The present paper is distinguished from its predecessors mainly in the specificity of 
the predictions (representative 99% C.L. predictions for each of the six currently acceptable 
neutrino oscillation scenarios) and in the much larger number of measurable quantities for 
which we now make predictions. 

Recent review articles summarize clearly the present state of neutrino physics |T3|JT^] and 



neutrino oscillation experiments and theory p,|15|dl7l1- Three and four flavor solar neutrino 
oscillations are discussed in Refs. [fLlj| , |l9|l and references cited therein. The fundamental 
papers upon which all of the subsequent solar neutrino oscillation work is based are the 



initial study of vacuum oscillations by Gribov and Pontecorvo |20[ and the initial studies 
of matter oscillations (MSW) by Mikheyev, Smirnov, and Wolfenstein [J2T|. In addition 
to the by-now conventional scenarios of oscillations into active neutrinos, we also consider 
oscillations into sterile neutrinos P, |H||!Z2| -pS|] . 
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II. ESTIMATION OF UNCERTAINTIES 



In this section, we describe how we calculate the uncertainties for different predicted 
quantities. Since the interpretation of future experimental results depends upon the assigned 
uncertainties, we present here a full description of how we determine the errors that we use 
in the remainder of the paper. 

Let X represent the predicted quantity of interest, which may be, for example, the first 
or second moment of the recoil energy spectrum, the neutrino-electron scattering rate, the 
double ratio of neutral current to charged current rate, the double ratio of neutrino-electron 
scattering to charged current rate, or the difference between the day rate and the night 
rate. The method that we adopt is the same in all cases. We evaluate X with two different 
assumptions about the size or behavior of a particular input parameter (experimental or 
theoretical). The different assumptions are chosen so as to represent a definite number of 
standard deviations from the expected best-estimate. The difference between the values of 
X calculated for the two assumptions determines the estimated uncertainty in X due to the 
quantity varied. 

To clarify what we are doing, we illustrate the procedure with specific examples. We 
begin by describing in Sec. [11 A| how we calculate theoretical uncertainties and then we 
discuss the detector-related uncertainties in Sec. Ill JJ. 



A. Theoretical uncertainties 

We discuss in this subsection the uncertainties related to the 8 B neutrino energy spec- 
trum, the neutrino interaction cross sections, and the hep solar neutrino flux. The standard 
shape of the 8 B neutrino spectrum has been determined from the best-available experimental 
and theoretical information |27| . 

Figure |3] shows the recoil electron energy spectra calculated for neutrino-electron scat- 
tering and for charged current (absorption) on deuterium that were calculated using the 
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FIG. 3. The calculated standard recoil electron energy spectra in SNO. The figure shows the pre- 
dicted shapes, A(T), of the normalized recoil electron energy spectra that are obtained by assuming 
that no oscillations occur and by using the standard (undistorted) 8 B neutrino energy spectrum. 
The spectra are given as a function of the true electron kinetic energy in MeV, T ej true; not the 
apparent energy measured by the detector. The spectra shown do not include instrumental effects 
such as the finite energy resolution of the detector or uncertainties in the absolute energy scale. 
The dotted curve represents the recoil electron spectrum due to neutrino-electron scattering and 
the solid curve represents the electron spectrum produced by neutrino absorption (CC reactions) 
on deuterium. 

undistorted standard 8 B neutrino energy spectrum. The recoil energy spectra produced by 
neutrino-electron scattering and by neutrino absorption are very different. One can easily 
see from Fig. |3] how the location of the threshold for CC events at 5 MeV (before the peak) 
or at 8 MeV (after the peak) could give rise to different sensitivities to uncertainties in, 
e.g., the energy resolution function. This is one of the reasons why we have calculated in 
the following sections predicted values and uncertainties for two different thresholds. For 
neutrino-electron scattering, the energy distribution decreases monotonically from low to 
high energies; this uniform behavior decreases the sensitivity, relative to the absorption 
process, to some uncertainties. 
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TABLE II. Neutrino Interaction Cross Sections on 
Deuterium. The table gives, in units of 10~ 42 cm 2 , the neu- 
trino charged current (CC) and neutral current (CC) cross 
sections for deuterium. The cross sections calculated by dif- 
ferent authors (shown in column 1) have been averaged over 
a standard 8 B neutrino energy spectrum. For the CC reac- 
tions, a 5 MeV threshold for the recoil electron energy was 
assumed and the energy resolution function for SNO was ap- 



proximated by Eq. (10) and Eq. (11). The last column gives 
the calculated ratios of the NC to CC ratios. 



Authors 


CC 


NC 


NC/CC 


KN a 


0.979 


0.478 


0.488 


YHH b 


0.923 


0.449 


0.486 


EBL C 


0.889 







a Ref. [28|. 



b Ref. (29]. 



Ref. |(J and Ref. |L1| 



Two extreme deviations from the shape of the standard spectrum were also determined 
using the best-available information | Z7| ; these extreme shapes represent the total effective 
±3cr deviations. We calculate the quantities X that will be measured by SNO using the 
standard 8 B spectrum and the effective 3a different spectra and determine from the following 
formula the associated uncertainty due to the shape of the 8 B spectrum. Thus 



ax( B spectrum) = 6 [ | X(+3a spectrum) — X(—3a spectrum) 



(6) 



Table [TT] lists the three relatively recent calculations for the charged current absorption 
cross sections on deuterium, by Ying, Haxton, and Henley (YHH) El, by Kubodera and 
Nozawa (KN) [ESJ, and by Bahcall and Lisi (BL) [TT|; the YHH and KN calculations use 
potential models and BL used an effective range treatment. For the neutral current cross 
sections, only the YHH and KN cross sections are available. If the quantity X involves the 
neutral current, then we define the la uncertainty by evaluating 

<Tx(NC cross section) = | X(YHH cross section) — X(KN cross section) |. (7) 

Using the values given in Table [H|, we define analogous la uncertainties for quantities asso- 
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dated with the CC and the double ratio, [CC]/[NC]. 

There is no principle of physics that enables one to set a rigorous error estimate based 
upon the cross section calculations summarized in Table [TT]. As a practical and plausible 
estimate for this paper, we have used the average of the detailed Kubodera and Nozawa and 
Ying, Haxton, and Henley calculations as our best estimate and taken the difference between 
these two cross sections to be an effective la uncertainty (see also the discussion by Butler 
and Chen in Ref. [31 1). Experimental measurements with reactor anti- neutrinos are not 



yet sufficiently accurate to refine and choose between different theoretical calculations (see 
results in Ref. Had we adopted the Ellis, Bahcall, and Lisi effective range calculation 



as the lower limit instead of the Ying et al. result, we would have obtained for the CC 
an uncertainty of 9.7% instead of 5.8%. Earlier, Bahcall and Kubodera |33| estimated 



an effective 3a uncertainty of ±10% for the neutral current cross section by calculating 
cross sections with and without meson-exchange corrections, using different sets of coupling 
constants, and two different nuclear potentials. 

The nuclear fusion reaction that produces hep neutrinos cannot be calculated or measured 
reliably P^ , |5B| . The shape of the electron recoil energy distribution measured by Super- 



Kamiokande can be significantly influenced by the rare high-energy hep neutrinos PJ3~4}-|3T)|1 . 
In this paper, we need to evaluate the uncertainty in a variety of quantities X due to the 
unknown hep flux. We use the results given in the last column of Table 3 of Ref. [|J, 
which lists the range of hep fluxes that correspond to different oscillation solutions that lie 
within the 99% (~ 2.5a) C.L. allowed range. Given the range of listed hep fluxes, we make 
the plausible but not rigorous estimate that the effective la uncertainty in the hep flux 
is currently between and 20 times the nominal standard estimate of 2.15 x 10 3 cm _2 s _1 
(0.0004 the best-estimate 8 B flux). Therefore, we evaluate the uncertainty due to the increase 
of the hep flux above the nominal standard value from the following relation 

a x hep flux = | X(20 x <j>(kep, BP98)) - X(0 x <j>(hep, BP98)) |. (8) 

The uncertainty in the hep flux is asymmetric (negative fluxes are not physical). We calculate 
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the lower error by replacing 20 x <p(hep, BP98) in Eq. (^) by 1 x <f>(hep, BP98). The lower 
error corresponds to decreasing the hep flux to zero. The uncertainty in the hep flux does 
not dominate the error budget for any of the quantities we discuss. If the reader wishes to 
treat differently the hep flux uncertainty, this can be done easily by using the individual 
uncertainties in Table |m[ 

For the standard solar model (SSM), the nominal ratio of the hep neutrino flux to the 
8 B neutrino flux is 4 x 10~ 4 [||]. Of all the quantities we consider in this paper, the first and 
second moments of the electron recoil energy spectrum, which are discussed in Sec. [TH[ are 
most sensitive to the hep flux. For a nominal SSM hep flux, the first moment is shifted by 
3 x 10~ 4 relative to the first moment computed with a zero hep flux. The corresponding 
change for the standard deviation of the recoil energy spectrum is 2 x 10~ 3 . Thus the hep flux 
of the standard solar model is of negligible importance for all of the quantities we calculate 
in this paper. The hep neutrino flux will have a significant effect on the quantities computed 
here only if the flux exceeds the nominal standard value by at least an order of magnitude. 

Super-Kamiokande and SNO will obtain somewhat tighter constraints on the hep flux. 
Measurements of the seasonal variations of the 7 Be flux will test vacuum neutrino scenarios 
that have a small hep flux but an appreciable distortion of the Super-Kamiokande recoil 
energy spectrum [p7|l . Since <Jx{hep flux) is linearly proportional to the allowed range of the 
hep flux, a reduction in the allowed range by, for example, a factor of two will reduce the 
estimated value of ax{hep flux) by a factor of two. 

For neutrino-electron scattering, the situation is very different. The interaction cross 
sections can be calculated precisely including even the small contributions from radiative 



corrections. We use in this work the cross sections calculated in Ref. [pq] ; the uncertainties 
in these radiative corrections are negligible for our purposes. 

In the following sections, we often quote fractional uncertainties in percent. We define 
the fractional uncertainty to be the one sigma difference divided by the average of the two 
values used to obtain the error estimate. Thus the fractional uncertainty due to an increase 
in the poorly known hep flux is 
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2&x(hep flux) 

5X(he P flux) = 100 x |x(2Ox0( ^ Bpg8)) + X(Ox0 (/^, BP98))|- (9) 



B. Detector-related uncertainties 

There are important detector-related uncertainties that can only be determined by de- 
tailed measurements with the SNO detector and by careful Monte Carlo simulations. Per- 
haps the most dangerous of these uncertainties are the misidentification uncertainties, the 
incorrect classification of CC, ES, and NC events. These errors do not cancel in the double 
ratios discussed later in this paper, such as [ES]/[CC] and [NC]/[CC]. 

SNO is a unique detector. No other detector has previously separated the CC and the NC 
reactions. The only reliable way of estimating the effects of the confusion between different 
neutrino reactions is to use the full-scale Monte Carlo simulation that is under development 
by the SNO collaboration. Since the SNO collaboration will measure the NC rate in different 
ways, there will also ultimately be internal cross checks that will limit the error due to the 
NC contamination of the CC and the ES rates. The ES contamination of CC quantities like 
the spectrum distortion or the day-night effect is likely to be small, since neutrino-electron 
scattering is strongly peaked in the forward direction and is estimated to be detected at 
only ~ 0.1 the CC rate. Hopefully, misclassification errors will have only minor effects and 
will be well described by the SNO Monte Carlo simulations. But, the reader should keep in 
mind that the errors estimated in this paper are lower limits; they represent errors that we 
can estimate quantitatively without a large Monte Carlo simulation. We will not consider 
errors due to misclassification of neutrino event in the reminder of this paper. 

One can make reasonable guesses for other important experimental uncertainties using 
the experience gained from previous water Cherenkov solar neutrino experiments and pre- 
liminary Monte Carlo studies of how the SNO detector will perform. The most important 
of these quantities that need to be determined, together with their uncertainties, are the 
energy resolution, the absolute energy scale, the detector efficiencies (for energetic electrons 
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and for neutral current reactions), and the energy threshold for detecting CC events. In 
what follows, we will adopt the preliminary characterizations for these detector-related un- 
certainties used by Bahcall and Lisi [f[T]j . We now summarize briefly our specific assumptions 
for these uncertainties. 

Let T' e be the true electron recoil kinetic energy and T e be the kinetic energy measured 
by SNO. We adopt the resolution function R(T' e , T e ), 



a(Tl)y/2n 



exp 



2a(T>y 



(10) 



with an energy-dependent one-sigma width cr(Tg) given approximately by 

a(T' e ) = (1.1 ±0.11) J — T { „ MeV. 



10 MeV 



(11) 



We adopt a conservative estimate for the la absolute energy error of ±100 keV. We will 
assume, for illustrative purposes, that the threshold for detecting recoil electrons is a total 
energy of 5 MeV or 8 MeV. 

For specificity, we assume jTT| that the neutral current detection efficiency is 0.50 ±0.01 
and that the detection efficiency for recoil electrons above threshold is approximately 100%. 



C. Summary of uncertainties 

In this subsection, we present a convenient table that summarizes the estimated un- 
certainties for the different physical quantities that are discussed in detail in the following 
sections of the paper. It may be useful to the reader to refer back to this summary table 
from time-to-time while considering the detailed presentations. 

Table |T| shows the fractional uncertainties in percent that we have estimated for different 
measurable quantities. The quantities in the Table are defined in the following sections. The 
counting uncertainties are determined assuming that a total of 5000 events are measured 
in the CC mode; the number of NC and neutrino-electron scattering events are then about 
1219 and 458, respectively. For an 8 MeV electron energy threshold rather than the 5 MeV 
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TABLE III. Fractional uncertainties in percent for some quantities that are measur- 
able with the SNO detector. Here SX = lOOax / X. The different quantities are denned in the 
following sections: T and a are the first and second moments, respectively; [CC], [NC], and [ES] are 
the reduced charged current, neutral current, and neutrino-electron scattering rates, respectively; 
and [NC]/[CC] and [ES]/[CC] are the neutral current to charged current and neutrino-electron 
scattering to charged current double ratios, respectively. For CC reactions, a 5 MeV threshold was 
assumed for the energy of the recoil electrons. The statistical uncertainties are computed assuming 
5000 CC events, 1219 NC events, and 458 ES events. We assumed a NC detection uncertainty of 
la = 2%. We do not include uncertainties due to misclassification of neutrino events. 



Source 


ST 


5a 


8 [CC] 


S [NC] 


8 [ES] 


S [NC]/[CC] 


S [ES]/[CC] 


Energy resolution 


0.3 


1.4 


0.4 


~ 


0.1 


0.4 


0.3 


Energy scale 


0.8 


1.1 


1.5 


~ 


0.5 


1.5 


1.0 


B8 spectrum 


0.4 


0.8 


1.9 


1.6 


1.4 


0.3 


0.6 


Cross-section 


0.03 


0.15 


5.8 


6.4 


~ 


0.4 


5.8 


Statistics* 1 


0.35 


1.1 


1.4 


2.9 


4.7 


3.2 


4.9 


hep° 


0.8 


0.8 


2.3 


2.2 


1.6 


0.1 


0.7 


Total 


1.3 


2.4 


6.7 


7.4 


5.2 


3.6 


7.6 



a Not including background from other sources. 
b One a hep upper error. 



threshold used in computing Table p| , the statistical uncertainties would be increased by 
about 50% for the purely CC quantities like 5T, 5a, and [CC]. For quantities related to the 
ES rate, the statistical uncertainties would be increased by a factor of about 1.9 by raising 
the electron recoil energy threshold to 8 MeV. The statistical error is expected, for an 8 
MeV threshold, to dominate the uncertainty in the ES rate. 

There will be additional contributions to the statistical errors from background sources; 
these uncertainties can only be determined in the future from the detailed operational charac- 
teristics of the SNO detector. For example, the background from the CC events will increase 
the estimated statistical error for the neutrino-electron scattering events; the amount of the 
increase will depend upon the angular width of the peak in the v-e scattering function. We 
have not estimated uncertainties for the day-night asymmetry, A, denned by Eq. (|28|) since 
a detailed knowledge of the detector is required to estimate the small uncertainties in A. 



The errors due to the uncertainties in the hep flux are asymmetric. We show in Table QI 



only the upper limit uncertainties for hep. The lower limit uncertainties are negligibly small 
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for hep, since the standard model flux ratio for hep to 8 B is 0.0004. 

The actual background rates in the SNO detector are not yet known and may differ 
considerably from the rates that were estimated prior to the building of the observatory. We 
have therefore not attempted to include background uncertainties, although these may well 
be important for some of the quantities we calculate. 

For both the CC ratio of measured to standard model rate, [CC], and the similarly 



defined neutral current ratio, [NC], Table |TTl| shows that the absolute value of the neutrino 



cross section is the dominant source of uncertainty. This uncertainty almost entirely cancels 
out in the double ratio of ratios, [NC]/[CC]. The absolute energy scale and the value of the 
hep neutrino energy flux are the largest estimated uncertainties for the first moment of the 
CC recoil energy spectrum, (T). Counting statistics, assuming a total of 5000 CC events, is 
estimated to be the most important uncertainty for the neutrino-electron scattering ratio, 
[ES], and the neutral current to charged current double ratio [NC]/[CC]. 

In the subsequent discussion, we follow the frequently adopted practice of combining 
quadratically the estimated ex's from different sources, including theoretical errors on cross 
sections and on the hep flux. If the reader prefers to estimate the total uncertainty using a 
different prescription, this can easily be done using the individual uncertainties we present. 

III. THE SHAPE OF THE CC ELECTRON RECOIL ENERGY SPECTRUM 

In this section, we make use of the fact that solar influences on the shape of the 8 B 
neutrino energy spectrum are only of order 1 part in 10 5 |39[], i.e., are completely negligible. 
Therefore, we compare all of the neutrino oscillation predictions to the calculated results 
obtained using an undistorted neutrino spectrum inferred from laboratory data |27j . 



Figure shows as a solid line the calculated CC electron recoil energy spectrum that 
would be produced by an undistorted 8 B neutrino energy spectrum. The result shown in 
Fig. |] does not include instrumental effects such as the energy response of the detector, but 
best-estimates of the instrumental effects (see discussion in Sec. |11B[ ) are included in the 
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results given here and in the following sections. 

It is useful in thinking about the shapes of the different electron recoil energy spectra to 
consider the ratio, R(E e ), of the electron energy spectrum produced by a distorted neutrino 
spectrum to the spectrum that is calculated assuming a standard model neutrino energy 
spectrum pp|] . We define 



P(F n _ /( 8B )^(^) + f(hep)N hep (E e ) 

^) N fM( Ee ) + N^(E e ) [LZ) 

where f{ 8 B) is the ratio of the true 8 B neutrino flux that is created in the sun to the 
standard solar model 8 B neutrino flux, i.e., f( 8 B) = (f>( 8 B) true /(j)( 8 B) SSM . The quantity 
f(hep) is similarly defined as the ratio of true to standard solar model hep flux. Ng SM (E e ) 
is the number of events in a 0.5 MeV energy bin centered at E e and calculated for the SSM 
8 B neutrino flux without oscillations. Ns{E e ) is the same quantity with oscillations taken 
into account. Nfj^jf-(E e ) and Nh ep {E e ) are the corresponding numbers for the hep flux. We 
have included the instrumental effects as described in Sec. [11 B| . 

We determine the best-fit value of /( 8 B) and f(hep) for each pair of values of the oscil- 
lation parameters, Am 2 and sin 2 26, by comparing the theoretical predictions with the total 
rate and the recoil electron energy spectrum of Super-Kamiokande 

f B = f B (Am 2 ,sm 2 26), f hep = f hep (Am 2 , sin 2 26). (13) 

The uncertainties in the values of /( 8 B) and f(hep) are reflected in the allowed range of 
Am 2 and sin 2 26, but are not included explicitly in Table |TTT. 



Figure f| shows the ratio R(E e ) calculated for the six best-fit oscillation solutions. The 
values of R(E e ) are given in 0.5 MeV bins except for the last energy bin, where we include 
all CC events that produce recoil electrons with observed energies above 14 MeV. Only a 
few events (less than 1% of the total number of CC events) are predicted Q to lie above 
14 MeV since the 8 B neutrino energy spectrum barely extends beyond 14 MeV and the hep 
neutrinos, which extend up to 18.8 MeV, are expected to be very rare. 

Ultimately, SNO will measure the detailed shape of the CC recoil energy spectrum and 
compare the measurements with the full predictions of different oscillation scenarios, as 
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FIG. 4. The relative CC recoil energy spectra for six oscillation solutions. The figure shows the 
rates, R(E e ), predicted by different best-fit oscillation solutions as a function of the electron recoil 
energy, E e , divided by the rate predicted by the standard solar model and no neutrino oscillations. 
The oscillation solutions are described in Table |. The highest energy bin represents the average 
value of R(E e ) for electron energies between 14 MeV and 20 MeV. We use the best-fit hep flux 
found for each neutrino oscillation scenario, which affects the result in the highest energy bin. The 
range of hep fluxes is given in Table 3 of Ref. 0] . 

illustrated in Fig. |j. Since the neutrino oscillation parameters are continuous variables, 
there are in principle an infinite number of possible shapes to consider. However, much or 
most of the quantitative information can be summarized conveniently in the first and second 
moments of the recoil energy spectrum [12] and we therefore concentrate here on the lowest 
order moments. 

Throughout this section, we use the notation of Bahcall, Krastev, and Lisi [12|] (hereafter 
BKL97), who have defined the first and second moments (average and variance) of the 
electron recoil energy spectrum from CC interactions in SNO. The explicit expressions are 
given in Eqs. (11)— (17) of BKL97; they include the energy resolution function of the detector 
[see Eq. fllCf ) of this paper]. Unlike BKL97, we use as our default recoil energy spectrum 5 
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MeV total electron energy, rather than 5 MeV electron kinetic energy. (We also calculate 
the moments for an 8 MeV total electron recoil energy.) When we calculate for the same 
threshold as BKL97, our results for the no-oscillation solution agree to about 1 part in 10 4 . 
We use a threshold specified in terms of total electron energy because this variable has 
become the standard for experimentalists to specify their energy threshold. 

We denote by a subscript of "0" the standard value of quantities computed assuming 
no oscillations occur. In order to compare with the theoretical moments given here, the 
observed moments should be corrected for any dependence of the detection efficiency upon 
energy that is determined experimentally. 

If there are no oscillations, the first moment of the CC electron recoil kinetic energy 
spectrum is, for a 5 MeV total electron energy threshold: 

(T) = 7.422 x (1 ± 0.013) MeV, (14) 

where the estimated uncertainties (±96 keV) have been taken from Table III]. The result 
given in Eq. ( |I4"D applies for a pure 8 B neutrino spectrum. If one includes a hep neutrino 
flux equal to the nominal standard solar model value H , then the first moment is increased 
by 2 keV to 7.424 MeV. For an 8 MeV energy threshold, (T) = 9.117 for a pure 8 B neutrino 
energy spectrum and is increased by 3 keV by adding a nominal hep flux. 

The largest estimated contributions to the quoted error in Eq. (|14D arise from uncer- 
tainties in the energy scale and from the hep reaction, with smaller contributions from the 
width of the energy resolution function and the shape of the 8 B neutrino energy spectrum. 
The total error of the measured value is the same, within practical accuracy, whether or not 
one includes the statistical uncertainty for 5000 events. 

The first moment has the smallest estimated total error of all the quantities tabulated 
in Table [TJ. 

Table [TV] presents the best-estimates and the total range of the predictions for the six 
different two-flavor neutrino scenarios that are globally consistent with all of the available 
neutrino data. Figure 1 of Ref. [£| shows, at 99% CL, the allowed ranges of the neutrino 
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TABLE IV. The first moment, (T), of the electron recoil kinetic energy spectrum 
from CC interactions. If there are no oscillations, the expected value of the first moment is 
(T ) = 7.422 MeV for a 5 MeV total electron energy threshold and (T ) = 9.117 MeV for an 8 MeV 
energy threshold. The table shows for different neutrino oscillation scenarios the absolute shift, 
AT = (T) — (T)o, in keV of the first moment of the electron recoil energy spectrum. Results are 
given for both a 5 MeV and an 8 MeV threshold energy. The different columns give the best-fit 
shift as well as the minimum and maximum shifts at 99% C.L. 



Scenario 


(AT) b . f . 


(AT) min 


(AT) max 


(AT) b . f . 


(AT) min 


(AT) max 




keV 


keV 


keV 


keV 


keV 


keV 




5 MeV 


5 MeV 


5 MeV 


8 MeV 


8 MeV 


8 MeV 



LMA 8 -115 34 4 -35 15 

SMA 218 50 341 66 15 105 

LOW 12 -17 63 7 -5 25 

VAC S 283 -80 576 122 40 227 

VAC L 21 -152 214 236 -54 358 

Sterile 164 41 265 51 13 83 



oscillation parameters of the first five neutrino scenarios listed in Table [IV[ The abbrevia- 
tions LMA, SMA, and LOW represent three MSW solution islands and the abbreviations 
VACs and VACl represent the small-mass and large mass vacuum oscillation solutions, all 
for oscillations into active neutrinos. The MSW Sterile solution has values for the mixing 
angle and the square of the mass difference that are similar to the active SMA solution (see 
discussion in Ref. (4fj). 

For a 5 MeV electron energy threshold, the predicted shifts in the first moment, AT = 
(T) — (T) , range from —152 keV to +576 keV. The calculational uncertainties and the 
measurement uncertainties estimated from the expected behavior of SNO, ±96 keV, are 
considerably smaller than the total range of shifts, 711 keV, predicted by the currently 
allowed set of oscillation solutions. The shift in the first moment may be measurable if 
either the SMA, VAC S , VAC L , or MSW Sterile solutions are correct. For the LMA and 
LOW solutions, the predicted shifts in the first moment may be too small to obtain a very 
significant measurement. 

A measurement of the first moment with an energy threshold of 5 MeV and a la accuracy 
in (T) of 100 keV or better will significantly reduce the allowed range of neutrino oscillation 
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solutions. Table [IV| shows that a measurement of (T) with an energy threshold of 8 MeV 
will be valuable, although it will provide a less stringent constraint than a measurement 
with a lower threshold. For an 8 MeV threshold, the currently allowed range is only 412 
keV, almost a factor of two less than the range currently allowed for a 5 MeV threshold. 
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FIG. 5. The fractional shift in the average electron recoil energy. The figure shows 5T [defined 
in Eq. fllq)], the fractional change in the average electron recoil energy, (T), for the six currently 
allowed neutrino oscillation solutions. The solid error bars represent the 99% C.L. for the allowed 
regions of the six currently favored neutrino oscillation solutions 0] . For an undistorted 8 B neutrino 
energy spectrum, the average recoil energy is denoted by To. The dashed error bar labeled "Measure 
3cr" represents the uncertainty in interpreting the measurements according to the estimates in 



Table III, which include the energy resolution, energy scale, B neutrino energy spectrum, neutrino 
cross section, counting statistics, and the hep flux. The results are calculated assuming a 5 MeV 
threshold for the CC reaction. 



Figure [| shows, for a 5 MeV electron energy threshold, the range of the fractional shift 
in percent of the first moment, 



5T = AT/(T) , 



(15) 



for all six of the oscillation solutions. The results are compared with the no-oscillation 
solution, 5T = 0. The estimated experimental uncertainty in 5T is about 1% (see Table [TTT]). 
Only the VACs solutions predict, for about half of their currently allowed solution space, a 
deviation of (T) from the no-oscillation value by more than 3a. The MSW sterile solution 
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TABLE V. The standard deviation, (er 2 ) 1 / 2 , of the electron recoil energy spectrum. 

If there are no oscillations, the expected value of the standard deviation of the electron recoil 
energy spectrum is <7o = (c 2 )c/ 2 = 1.852 MeV for a 5 MeV total electron recoil energy threshold 
and 1.240 MeV for an 8 MeV energy threshold. The table shows for different neutrino oscillation 
scenarios the absolute shift of the standard deviation, Aa = (cr 2 ) 1 / 2 — (a 2 )^ 2 , in keV for both a 5 
MeV total electron recoil energy threshold and an 8 MeV energy threshold. 



Scenario 


(A*) b . f . 


( Af7 )min 


(Aa) max 


(Aff)bj. 


(Acr) min 


(Acr) max 




keV 


keV 


keV 


keV 


keV 


keV 




5 MeV 


5 MeV 


5 MeV 


8 MeV 


8 MeV 


8 MeV 



LMA 3 -19 9 4 -13 11 

SMA 23 1 38 23 3 34 

LOW 3 -3 13 2 -3 10 

VAC S 70 11 136 44 9 76 

VAC L 127 -29 199 160 -40 212 

Sterile 19 2 32 14 -4 36 



predicts a shift in the first moment that is at most 2.7a from the no-oscillation case; this 
seems like a small shift, but it is notoriously difficult to identify measurable indications of 
sterile neutrinos that are different from a reduction in the total 8 B solar neutrino flux |26l . 



Table |V| presents the predicted shifts in the standard deviation of the CC electron re- 
coil energy distribution (i.e., the square root of the second moment). The calculated no- 
oscillation value is 

a = ( a 2 )l /2 = 1.852(1 ± 0.024) MeV, (16) 

for a 5 MeV total electron recoil energy threshold and o"o = 1.240 MeV for an 8 MeV 
threshold. The estimated uncertainties in Eq. ( |16D are taken from Table |IT| and Table III 



of Ref. [ I2| . The result given in Eq. fllB]) is for a pure 8 B neutrino spectrum. If a hep flux 



equal to the standard solar model value is included, the value of <7o is increased by 4 keV 
to 1.856 keV. For an 8 MeV threshold, o"o is increased by 6 keV to 1.246 MeV by adding a 
nominal standard hep flux. 

Shifts in the standard deviation caused by neutrino oscillations will be difficult to measure 
since the spread in the predicted shifts for a 5 MeV threshold is only from —29 keV to +199 
keV, while the estimated calculational and non-statistical measurement uncertainties are 
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±91 keV. Thus the total range of the predicted shifts is less than three standard deviations 
of the estimated non-statistical uncertainties. For most of the oscillation scenarios, the shift 
in u predicted for an 8 MeV threshold is even smaller than for a 5 MeV threshold. 

It will be useful to measure the standard deviation of the recoil energy spectrum in order 
to test the prediction that the observed value will be close to the undistorted value of <t . 



IV. THE CHARGED CURRENT RATE 

In this section, we summarize the results from Ref. [[| on the expected range of predictions 
for the charged current (neutrino absorption) rate [see Eq. ([!])]. 

In accordance with Eq. @, we define the reduced CC neutrino- absorption ratio [CC] by 

r , _ (Observed CC Rate) 

[ J = (Standard CC Rate) ' ^ ' 

If the standard solar model is correct and there are no neutrino oscillations or other non- 
standard physics processes, then 



[CC] = 1.0 
= 1.0 



1.0 ±0.058 a ±0.019 b ±0.004 c ± 0.015 d to^o 

i n +0.067" 
l - v -0.061 



(18) 



The uncertainties in the standard solar model flux || do not affect [CC] since we fix the 
absolute flux for each set of neutrino parameters by fitting to the Super-Kamiokande total 
rate and recoil energy spectrum (see discussion following Eq. |13|). 

The reduced v-e scattering rate has been measured by Super-Kamiokande to be about 
0.475 ± 0.015 for a threshold of 6.5 MeV and the reduced v-e scattering ratio is predicted 
to be approximately the same for the expected SNO energy thresholds (see Table |V11| and 



Ref. \WT 



The non-statistical uncertainties shown in Eq. (|T%| ) result from (cf. Table |111| ): (a) the 
difference between the Ying, Haxton, and Henley [PJ and Kubodera-Nozawa cross sections 
neutrino cross sections, (b) the shape of the 8 B neutrino energy spectrum E7[, (c) the 
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TABLE VI. The Charged Current Absorption Ratio, [CC]. The table presents, [CC], 
the ratio of the observed neutrino absorption rate on deuterium to the standard model absorption 
rate, cf. Eq. (|l7|). The results are tabulated for different neutrino oscillation scenarios and for 
two different thresholds of the total electron recoil energy, 5 MeV (columns two through four) and 
8 MeV (columns five through seven). The second (fifth) column gives the best-fit value, and the 
third (sixth) and four (seventh) columns give the minimum and maximum values. 



Scenario 


[CC] 


[CC] 


[CC] 


[CC] 


[CC] 


[CC] 




b.f. 


max 


min 


b.f. 


max 


min 




5 MeV 


5 MeV 


5 MeV 


8 MeV 


8 MeV 


8 MeV 



LMA 0.35 0.40 0.29 0.35 0.40 0.29 

SMA 0.39 0.46 0.32 0.43 0.47 0.38 

LOW 0.38 0.40 0.35 0.38 0.40 0.35 

VAC S 0.39 0.45 0.31 0.44 0.50 0.34 

VAC L 0.42 0.44 0.40 0.40 0.47 0.35 

Sterile 0.48 0.49 0.47 0.50 0.55 0.48 



energy resolution function, (d) the absolute energy scale, and, the last term, the uncertain 
hep neutrino flux. The total uncertainty in the charged current ratio [CC] is dominated by 
the uncertainty in the CC absorption cross section. 

The most important question concerning the CC rate in SNO is the following: Is the 
reduced CC rate less than the reduced neutrino-electron scattering rate? If the reduced CC 
rate is measured to be less than the reduced v-e scattering rate, then this will be evidence 
for neutral currents produced by v^ or v T which appear as a result of neutrino oscillations. 
The 3a uncertainty is about 20% above the expected value of 0.48 (based upon the Super- 
Kamiokande v-e scattering measurement). Inspecting for all six of the currently favored 



oscillation solutions the range predicted for the double ratio [CC] (as shown in Table VI or 
Fig. 2 of Ref. [[|), we estimate that there are very roughly equal odds that the measured 
value of [CC] will lie three or more a below the no-oscillation value. However, for the VACl 
and MSW Sterile solutions, the predicted range for [CC] does lie within 3a of the value 
expected on the basis of the no-oscillation hypothesis. 
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V. THE NEUTRAL CURRENT RATE 



We discuss in this section the expected range of predictions for the neutral current rate 
[see Eq. ©]. 

If the standard solar model is correct and if there are either no neutrino oscillations or 
oscillations only into active neutrinos, then the reduced NC neutrino-absorption rate, [NC], 
defined by 



— (Observed NC Rate) 
^ J = (Standard NC Rate) ' ^ ' 



will satisfy 



[NC] = 1.0 ±0.060 a ±0.016 b ±0.02 c _ ll000 _„ 1( , 

= 1.0 ±0.07 +°;} 8 . (20) 

The non-statistical uncertainties shown in Eq. (|20|) result from (cf. Table p|) : (a) the 
difference between the Ying, Haxton, and Henley and Kubodera-Nozawa cross sections 



neutrino cross sections, (b) the shape of the 8 B neutrino energy spectrum P7[ , and (c) 



the uncertainty in the neutral current detection efficiency The next to last term in Eq. 
represents the uncertainty in the hep neutrino flux. 

The last term in Eq. (p0| ) represents the uncertainty in the BP98 standard 8 B flux ||. 
In our method, this uncertainty only appears if there are no neutrino oscillations or other 
new physics. If there are neutrino oscillations, we determine the ratio, f( 8 B). of the best-fit 
neutrino flux to the standard model flux as described in Sec. p] following Eq. (fD and in 
Sec. |TJ following Eq. (0) and Eq. (|13]). 



The total uncertainty in determining experimentally the neutral current ratio [NC] is 
dominated by the uncertainty in the NC absorption cross section and the total uncertainty 
in interpreting the neutral current measurement is dominated by the uncertainty in the 
predicted solar model flux. 

If one assumes that only oscillations into active neutrinos occur, then it will be possible 
to use the neutral current measurement as a test of the solar model calculations. The cross 
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section uncertainty for the NC reaction is about 33% (38% ) of the upper (lower) estimated 
uncertainty in the solar model flux. 

If the SMA Sterile neutrino solution 0,^,0 is correct, then for the global solutions 
acceptable at 99% C.L., 

[NC] = 0.465 ± 0.01 . (21) 

Unfortunately, this result is within about 3<r of the result expected if there are oscillations 
into active neutrinos, when one includes the solar model uncertainty shown in Eq. (|20|) . 

VI. THE NEUTRINO-ELECTRON SCATTERING RATE 

In this section, we present the predictions for the electron-scattering rate, Eq. (|]), in 
SNO. The SNO event rate for this process is expected to be small, ~ 10% of the CC rate. 
Despite the relatively unfavorable statistical uncertainties, the measurement of the neutrino- 
electron scattering rate in SNO will be important for two reasons. First, the measurement of 
the electron scattering rate by SNO will provide an independent confirmation of the results 
from the Kamiokande and Super-Kamiokande experiments. Second, the neutrino-electron 
scattering rate can be combined with other quantities measured in SNO so as to decrease the 
systematic uncertainties and to help isolate the preferred neutrino oscillation parameters. 

If the best-estimate standard solar model neutrino flux || is correct and there are no 
new particle physics effects, then the reduced neutrino-electron scattering rate, [ES], will be 
measured to be 

l O LcLIlU.cLI u. v — e rvcLLcl 

where the non-statistical uncertainties are taken from Table |TJ. For a 6.5 MeV threshold and 
the experimental parameters of the Super-Kamiokande detector 0, [ES]sk = 0.475 ±0.015. 

The uncertainties summarized in Eq. ( |22] ) include all the uncertainties given in Table [Tl 1 
except for statistical errors. The uncertainties in the standard solar model flux M do not 
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TABLE VII. The Neutrino-Electron Scattering Ratio, [ES]. The table presents, [ES], 
the ratio of the observed neutrino-electron scattering rate to the standard model rate, cf. Eq. (|22|). 
The results are tabulated for different neutrino oscillation scenarios and for two different thresholds 
of the total electron recoil energy, 5 MeV (columns two through four) and 8 MeV (columns five 
through seven). The second (fifth) column gives the best-fit value, and the third (sixth) and four 
(seventh) columns give the minimum and maximum values. 



Scenario 


[ES] 


[ES] 


[ES] 


[ES] 


[ES] 


[ES] 




b.f. 


max 


min 


b.f. 


max 


min 




5 MeV 


5 MeV 


5 MeV 


8 MeV 


8 MeV 


8 MeV 



LMA 0.48 0.49 0.47 0.47 0.48 0.45 

SMA 0.47 0.48 0.46 0.50 0.51 0.48 

LOW 0.48 0.48 0.47 0.48 0.49 0.47 

VAC S 0.48 0.52 0.46 0.51 0.54 0.46 

VAC L 0.49 0.50 0.47 0.47 0.50 0.45 

Sterile 0.46 0.48 0.45 0.51 0.53 0.48 



affect [ES] since we fix the absolute flux for each set of neutrino parameters by fitting to the 
Super-Kamiokande total rate and recoil energy spectrum (see discussion following Eq. p~3| ) . 
The dominant non-statistical uncertainties are from the value of the hep flux and the shape 
of the 8 B neutrino energy spectrum. For the first five or ten years of the SNO operation, 
the overall dominant error in the determination of [ES] is expected to be statistical: > 5% 
after the first 5000 CC events. 



Table |VII| gives, for two different energy thresholds, the values of the reduced neutrino- 
electron scattering rate, [ES], that are predicted by the currently favored oscillation sce- 
narios [f| . Not surprisingly, the values of [ES] cluster around the ratio measured by the 
Super-Kamiokande experiment [ES]sk = 0.475 ± 0.015 0. The global constraints imposed 
by the different experiments result in some cases in the spread of the currently favored 
predictions being less than or of the order the total spread in the Super-Kamiokande rate 
measurement. 

Figure |6] compares the oscillation predictions for [ES]sno versus [ES]s upe rK and the no- 
oscillation solution. We only show the predictions for a 5 MeV threshold for the total electron 
energy since the results are similar for an 8 MeV threshold (see Table |V11| ). The solid error 
bars shown in Fig. reflect the range at 99% C.L. of the globally allowed solutions that are 
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FIG. 6. The reduced neutrino-electron scattering rate, [ES]. The figure shows the reduced neu- 
trino-electron scattering rate [cf. Eq. (|2^)] from the standard model value of 1.0 (no oscillations). 
The Super-Kamiokande result (see Ref. ||) is also shown. The solid error bars represent the 99% 
C.L. for the allowed regions of the six currently favored neutrino oscillation solutions 0. The 
dashed error bar labeled "Measure 3<r" represents the uncertainty in interpreting the measure- 



ments according to the estimates in Table III, which include the energy resolution, energy scale, 
8 B neutrino energy spectrum, neutrino cross section, counting statistics, and the hep flux. The os- 
cillation predictions are calculated assuming a 5 MeV threshold for the total electron recoil energy. 



fit to all the available neutrino data 



VII. THE NEUTRAL CURRENT TO CHARGED CURRENT DOUBLE RATIO 



In this section, we present predictions for the ratio of neutral current events (NC) to 
charged current events (CC) in SNO. The most convenient form in which to discuss this 
quantity is obtained by dividing the observed ratio by the ratio computed assumed the 
correctness of the standard electroweak model (SM). This double ratio is defined by the 



relation 10 



[NC] ((NC) bs/(NC) 
[CC] ((CCWCCC) 

The ratio [NC]/[CC] is equal to unity if nothing happens to the neutrinos after they are 
produced in the center of the sun (no oscillations occur). Also, [NC]/[CC] is independent of 
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TABLE VIII. Neutral Current to Charged Current Double Ratio. The table presents 
the double ratio, [NC]/[CC], that is defined by Eq. (|23|). The results are tabulated for different 
neutrino oscillation scenarios and for two different thresholds of the total electron recoil energy 
used in computing the CC ratio, 5 MeV (columns two through four) and 8 MeV (columns five 
through seven). 



Scenario 


[NC]/[CC] 


[NC]/[CC] 


[NC]/[CC] 


[NC]/[CC] 


[NC]/[CC] 


[NC]/[CC] 




b.f. 


max 


min 


b.f. 


max 


min 




5 MeV 


5 MeV 


5 MeV 


8 MeV 


8 MeV 


8 MeV 



LMA 3.37 5.15 2.27 3.36 5.13 2.30 

SMA 2.53 4.11 1.24 2.28 3.50 1.21 

LOW 2.71 3.39 2.30 2.69 3.35 2.26 

VAC S 2.67 4.63 1.37 2.33 3.73 1.27 

VAC L 1.90 2.15 1.53 2.01 2.45 1.53 

Sterile 0.96 0.99 0.94 0.88 0.97 0.82 



all solar model considerations provided that only one neutrino source, 8 B, contributes signif- 
icantly to the measured rates. Finally, the calcuiational uncertainties due to the interaction 
cross sections and to the shape of the 8 B neutrino energy spectrum are greatly reduced by 
forming the double ratio (see Table |m| ). 

Table |VIII| presents the calculated range of the double ratios for the oscillation solutions 
that are currently allowed at 99% CL 0. The table gives the best-fit values for [NC]/[CC] 
as well as the maximum and minimum allowed double ratios for a total electron energy 
threshold for the CC reaction of 5 MeV and separately for a CC threshold of 8 MeV. 

Figure |7| compares the predicted values of [NC]/[CC] with the no-oscillation value of 1.0. 
The results are shown for a 5 MeV CC threshold and for an 8 MeV CC threshold. The 
estimated [see Eq. (p4[)j non-statistical errors are smaller than the black dots indicating the 
best-fit points in Fig. [7[ 

The best-fit values for the double ratio for oscillations into active neutrinos range between 
1.9 and 3.4 (2.0 and 3.4) for a 5 MeV (8 MeV) CC threshold. The maximum predicted values 
for [NC]/[CC] exceed 5.1. For active neutrino oscillations, the minimum values for the double 
ratio are achieved by the SMA and the VACl solutions; they are 1.2 and 1.5, respectively. 

The sterile neutrino solutions predict a double ratio in a band that is separate from all 
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FIG. 7. The neutral current to charged current double ratio. The double ratio, [NC]/[CC] is 
defined by Eq. (^). The standard model value for [NC]/[CC] is 1.0. Figure [?|a shows, for a 5 
MeV threshold for the CC measurement, the predicted double ratio of Neutral Current to Charged 
Current for different neutrino scenarios. Figure 0b shows the same ratio but for an 8 MeV CC 
threshold. The solid error bars shown represent the 99% C.L. for the allowed regions of the six 
currently favored neutrino oscillation solutions 0. The dashed error bar labeled "Measure 3<t" 



represents the uncertainty in interpreting the measurements according to the estimates in Table 
which include the energy resolution, energy scale, 8 B neutrino energy spectrum, neutrino cross 
section, counting statistics, and the hep flux. 
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the active oscillation solutions, namely, 0.9 to 1.0. The physical reason that the double 
ratio for sterile neutrinos is less than 1.0 is that in the SMA solution (for active or sterile 
neutrinos) the probability that a solar v e survives as a v e decreases with energy (see, e.g., 
Fig. 9 of Ref. |J). In the sterile neutrino case, if the v e oscillates to another state it does 
not interact. Since the NC threshold is 2.2 MeV and the observational threshold for CC 
events is likely to be 5 MeV or above, the smaller survival probability at low energies more 
strongly affects the average NC rate than the average CC rate. 
The standard model value for [NC]/[CC] is 



[NC] 



1.0 ±0.004 a ±0.003 b ±0.004 c ± 0.015 d ± 0.02 c 



[CC] 

= 1.0 ±0.026. (24) 

The uncertainties, all non-statistical, shown in Eq. ( p4|) result from: (a) the difference be- 
tween the Ying, Haxton, and Henley |29] and Kubodera-Nozawa cross sections []28| , (b) the 



shape of the 8 B neutrino energy spectrum ^7j, (c) the energy resolution function, (d) the 
absolute energy scale, and (e) the NC detection efficiency Comparable, but small, contri- 
butions are made by the cross section uncertainties, the uncertainties in the shape of the 
neutrino energy spectrum, and the uncertainty in the energy resolution function. The ab- 
solute energy scale and the NC detection efficiency are expected to contribute even less to 
the errors. 

One of the principal uncertainties in interpreting the electron recoil energy spectrum 
is the poorly known value for the flux of the extremely rare hep neutrinos |35],|36|. The 
uncertainty in the hep flux can also affect the otherwise robust measurement of [NC]/[CC] 
(see Table [Hl|) . We have recalculated the value of [NC]/[CC] for a hep flux that is 20 times 
larger than the nominal standard model flux ||. We find [cf. Eq. (||) for the calculational 
prescription] 

[NC] / 4>(hep) \ 

[CC] " hep \20 (f>{hep, BP98)J' 1 ] 

where eh ep = 0.0005 for a 5 MeV threshold on the CC events and 0.017 for an 8 MeV 

CC threshold. For an 5 MeV CC threshold, there is an accidental cancellation of the 

31 



contributions to the neutral current ratio, [NC], and to the charged current ratio, [CC], 
so that the net value of e/i ep (5 MeV) is very small. But, for an 8 MeV threshold, the hep 
flux causes an uncertainty, 1.7%, that is larger than the combined contribution from all the 
other known uncertainties except possibly counting statistics [cf. Eq. ([M]) and Eq. and 
Table [ITJ . Of course, the gain at lower energies due to the reduction in the uncertainty from 
the hep neutrinos may be more than offset by the increased uncertainty due to background 
events. Fortunately, SNO is expected to be able to measure or to place strong limits on the 
hep flux within the first full year of operation 0]. 

We have not included the statistical uncertainties in the calculational error budget of 
Eq. fl2~4]). It seems likely that statistical errors will dominate over calculation errors, at 



least in the first several years of operation of SNO (see Table |III| ). The CC rate may be 
in the range of 3000 to 4000 events per year. The NC event rate in the detector will be 
about a factor of two smaller and the NC detection rate will be further decreased by the NC 
detection efficiency that may be of order 50%. Thus statistical errors in the NC rate, the 
uncertainty (~ 2%, see Table [TTTD in the NC detection rate, and the uncertainties in the hep 
flux, will probably be the limiting factors in determining the accuracy of the experimental 
measurement of [NC]/[CC]. 

The small calculational error [see Eq. (p4|)1 for [NC]/[CC], combined with the relatively 
large differences between the no-oscillation and the oscillation values for active neutrinos 
(Table |VIII|) , makes the double ratio an ideal 'smoking-gun' indicator of oscillations. 

VIII. THE ELECTRON-SCATTERING TO CC DOUBLE RATIO 

In this section, we present results for the double ratio of neutrino-electron scattering to 
CC events. This ratio is defined, by analogy with the NC to CC double ratio [see Eq. (p3])1, 
by the expression 

[ES] _ ((ES) bs/(ES) 



[CC] ((CC) bs/(CC) 



(26) 
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TABLE IX. Neutrino-Electron Scattering to Charged Current Double Ratio. The 

table presents the double ratio, [ES]/[CC], that is defined by Eq. (f26|). The results are tabulated 
for different neutrino oscillation scenarios and for two different thresholds of the total electron 
recoil energy for both the scattering and the CC reactions, 5 MeV (columns two through four) and 
8 MeV (columns five through seven). 



Scenario 


[ES/CC] 


[ES/CC] 


[ES/CC] 


[ES/CC] 


[ES/CC] 


[ES/CC] 




b.f. 


max 


min 


b.f. 


max 


min 




5 MeV 


5 MeV 


5 MeV 


8 MeV 


8 MeV 


8 MeV 



LMA 1.37 1.65 1.20 1.36 1.63 1.19 

SMA 1.20 1.43 1.03 1.29 1.57 1.05 

LOW 1.27 1.37 1.20 1.26 1.37 1.20 

VAC S 1.25 1.58 1.03 1.16 1.38 1.02 

VAC L 1.15 1.20 1.08 1.18 1.29 1.05 

Sterile 0.96 0.99 0.92 1.06 1.08 1.02 



The double ratio [ES]/[CC] has some of the same advantages as the double ratio 
[NC]/[CC], namely, independence of solar model considerations and partial cancellation 
of uncertainties. In fact, the [ES]/[CC] double ratio has the additional advantage that the 
same detection process is used for the recoil electrons from both the scattering and the CC 
reactions. For the [NC]/[CC] double ratio, different techniques are used to determine the 
two rates and this increases the systematic measurement uncertainty in the ratio. 

The standard model value for [ES]/[CC] is 

jS = 1.0 ±0.06 a ±0.006 b ±0.003 c ± 0.01 d 
[CCJ 

= 1.0 ±0.06, (27) 



where the non-statistical uncertainties shown in Eq. (27 ) result from: (a) the difference 



between the Ying, Haxton, and Henley j29] and Kubodera-Nozawa cross sections [28 



neu- 



trino cross sections, (b) the shape of the 8 B neutrino energy spectrum P7| , (c) the energy 
resolution function, and (d) the absolute energy scale. The upper limit hep uncertainty is 
small, 0.007 (see Table [IT]). 

The uncertainty in the double ratio [ES]/[CC] is dominated by the uncertainty in the 
CC absorption cross section and by statistical errors (~ 4.9% after 5000 CC events, see 
Table |ITJ) that are not included in Eq. fl27|) . 
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Table |IX| presents the calculated range of [ES]/[CC] for the oscillation solutions that are 
currently allowed at 99% CL J4]]. The table gives the best-fit values for [ES]/[CC] as well as 
the maximum and minimum allowed double ratios for a total electron energy threshold (for 
both reactions) of either 5 MeV or 8 MeV. The range of ratios predicted by oscillations into 
active neutrinos is 1.03 to 1.65, much smaller than the range (1.24 to 5.1) predicted for the 
[NC]/[CC] double ratio. 

Figure [| shows the values of [ES]/[CC] predicted by the different oscillation solutions. 
Comparing Fig. § and Fig. [7], one can see that the neutral current to charged current ratio 
is a more sensitive diagnostic of neutrino oscillations than is the electron scattering to 
charged current ratio. The difference from the no-oscillation solution is much greater for 
the [NC]/[CC] double ratio than it is for the [ES]/[CC] double ratio. In addition, there 
are expected to be many more detected NC events than neutrino-electron scattering events. 
Also, the cross section uncertainties largely cancel out of the ratio [NC]/[CC], whereas the 
uncertainty in the CC cross section is an important limitation in interpreting the ratio 
[ES]/[CC]. 



IX. THE DAY-NIGHT EFFECT 

We discuss in this section the difference between the event rate observed at night and 
the event rate observed during the day. For MSW solutions, the interactions with matter of 
the earth can change the flavor content of the solar neutrino beam and cause the nighttime 
and daytime rates to differ. This effect has been discussed and evaluated by many different 



authors, including those listed in Ref. ||42j| . 

We concentrate here on the difference, An-d, between the nighttime and the daytime 
rates, averaged over one year. The formal definition of An-d is 

[Night -Day] 

An ~ D " 2 [Night + Day]' (28) 
In what follows, we shall use An~d to refer to the charged current reaction. When we want to 
consider the quantity defined by Eq. (|28|) for the neutral current, we shall write An-d(NC). 
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FIG. 8. The neutrino-electron scattering to charged current double ratio. The double ratio, 
[ES]/[CC] is defined by Eq. (|26|). Figure [8|a shows, for a 5 MeV electron recoil energy threshold, 
the predicted double ratio of Neutrino-Electron Scattering to Charged Current for different neutrino 
scenarios. The solid error bars represent the 99% C.L. for the allowed regions of the six currently 
favored neutrino oscillation solutions 0. The dashed error bar labeled "Measure 3cr" represents 



the uncertainty in interpreting the measurements according to the estimates in Table 111, which 
include the energy resolution, energy scale, 8 B neutrino energy spectrum, neutrino cross section, 
counting statistics, and the hep flux. Figure [?]b shows the same ratio but for an 8 MeV threshold. 
Figure |8[ which has a vertical scale extending from 0.5 to 3.0, should be compared with Fig. 
which has a vertical scale extending from 0.5 to 7.0. The difference from the no-oscillation solution 
is much greater for the [NC]/[CC] double ratio than it is for the [ES]/[CC] double ratio. 
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We begin by discussing in Sec. [IX A| the apparent day-night effect that arises solely from 



the eccentricity of the earth's orbit and the inclination of the earth's axis (the existence of 
seasons), and then we discuss in Sec. |1X B| the day-night effect for the CC reaction and in 



Sec. [IX C| the day-night effect for the NC reaction due to oscillations. 

A. The No-Oscillation day-night effect 

In the absence of neutrino oscillations, there is a geometrical day-night effect that we 
have not seen discussed in previous publications. This No-Oscillation (NO) effect is caused 
by the ellipticity of the earth's orbit and by the fact that, in the northern hemisphere, nights 
are longer (days are shorter) in winter when the earth is closer to the sun. Thus the average 
over the year of the nighttime rate will be larger than the annual average of the daytime 
rate for all detectors located in the northern hemisphere. 

We find that the No-Oscillation (NO) day-night effect is 

A%°_ D = 0.0094 (SNO), (29) 

A^°_ D = 0.0066 (SK), (30) 

A^°_ D = 0.0082 (GranSasso), (31) 

^n-d = 0.0088 (Homestake), (32) 

for the locations of the SNO, Super-Kamiokande, Gran Sasso, and Homestake detectors. 

The No-Oscillation effect is purely geometrical; it is independent of neutrino energy and 
independent of neutrino flavor. The magnitude of the NO effect is the same for the CC, ES, 
and NC reactions. The numerical results given in Eq. fl30|)- Eq. fl32]) can also be obtained 
from the following easily- derived relation, which makes clear the seasonal aspect of the NO 
effect: 

jNo Osc. r) (^niax ~ ^min) , qq n 

where e = 0.0167 is the eccentricity of the earth's orbit and t max and t m i n are, respectively, 
the length of the longest and the shortest nights in the year at the location of the detector. 
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In what follows, we remove the No-Oscillation day-night effect before presenting the 
predictions of an additional day-night effect that is due to neutrino oscillations. More 
precisely, we calculate the day-night effect assuming that the neutrino flux from the sun is 
constant throughout the year. The effects that we discuss in Sec. |1X B| and in Sec. |1X C are 
due to neutrino mixing. 

Experimental results can easily be analyzed so as to remove the NO day-night effect. 
All that is required is to multiply the number of events in each time bin (At <C I year) 
by the ratio [r(t)/(I A.U.)] 2 , where r(t) is the average earth-sun distance in that time bin 
and 1 A.U. is the annual average earth-sun distance. This is the procedure adopted by the 
Super-Kamiokande collaboration [||]. 

Even after these corrections, there is a residual day-night effect for vacuum oscillations. 
In this case, the day-night effect is due to the variation of the survival probability as a 
function of the distance between the earth and the sun and the fact that in the northern 
hemisphere the longest nights occur when the earth is closest to the sun. We are not 
aware of any previous discussions of the day-night effect for vacuum oscillations. For MSW 
oscillations, the day-night effect is caused by neutrino flavor changes during propagation in 
the earth. 



B. The CC day- night effect 

Table [X] and Fig. || present the range of predicted percentage differences between the 
average rate at night and the average rate during the day [i.e., fOO x An-d of Eq. (§8|)]. 
The calculated predictions are given for a 5 MeV and an 8 MeV CC electron recoil energy 
threshold. 

For vacuum oscillations, the day-night effect is due to the dependence of the survival 
probability upon the earth-sun distance. The predicted day-night effect for vacuum oscilla- 
tions is small in all the cases shown in Table [X] and in Fig. |9|. 

For most of the MSW oscillation solutions, the predicted day-night differences are only 
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of order a few percent. However, for the LMA solution, the predicted difference can reach 
as high as 29% for a 5 MeV threshold (32% for an 8 MeV threshold). There are also rather 
large differences, in excess of 10%, that are possible for the SMA and LOW solutions. 

At first glance, one might think that such large day-night differences will be easy to 
measure. In fact, there are important systematic uncertainties that have to be taken into 
account in making sure that the relative sensitivities to the day and the night rates are prop- 
erly evaluated ||. Even the purely statistical uncertainties are very significant because the 
day-night difference, An-d, is the difference between two comparably sized large numbers. 
Thus the fractional statistical uncertainty after accumulating a large number, N, of counts 
at night (and a roughly equal number during the day) is 



The fact that can be a small number makes a multi-sigma statistical measurement 

of the day-night effect difficult. The careful analysis of the day-night effect by the Super- 



surement of A N _ D . Using more than 800 effective days of operation of the SuperK detector 
with total night time counts of N ~ 5900 (~ 11, 200 total events) the precision obtained by 
the Super-Kamiokande collaboration is A^k = 0.065 ±0.03, i.e., AA/A ~ 0.5. To accumu- 
late with SNO an equivalent number of CC events (~ 11,000 total events) may require of 
order three years or longer of operation. 

Why is the predicted effect in SNO (see also Ref. f43| ) so much larger than for Super- 
Kamiokande? The reason is that for neutrino-electron scattering the day-night effect is 
decreased relative to the pure CC mode by the contribution of the neutral currents. For 
the LMA solution, one can derive a simple quantitative relation between the CC day-night 
effect, A cc , and the ES day-night effect, A ES . Let the nighttime rate be proportional to 
Pn + r(l — -Pn), where Pn is the average (over energy) survival probability during the night 
and r is the average ratio of — e to v e — e scattering cross sections. Writing a similar 
expression for the daytime rate, it is easy to show that 




(34) 



Kamiokande collaboration || has demonstrated the practical difficulty of a precision mea- 
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TABLE X. The Predicted Night-Day Difference (in %) for the CC Rate. The table 
gives the percentage difference in the Night-Day CC rates, An_e>, defined by Eq. (p8|). For a 5 
MeV (8 MeV) threshold energy for recoil electrons, the second (fifth) column gives the best-fit 
value, and the third (sixth) and fourth (seventh) columns give the minimum and maximum values, 
allowed at 99% CL for different neutrino oscillation solutions (cf . Fig. 1 of Ref . 0] ) . 

Scenario b.f. min max b.f. min max 



LMA 12.4 +0.5 

SMA 1.6 -1.1 

LOW 4.7 +1.1 

VAC S 0.5 -0.1 

VAC L 0.3 -0.1 

MSW, Sterile -0.1 -0.4 



28.5 14.1 0.7 31.6 

12.3 1.2 -1.3 10.9 

13.5 4.0 0.8 11.8 

1.0 0.7 0.2 1.2 
0.5 1.0 -0.1 1.5 

1.1 -0.2 -0.7 1.1 



4 cc 



4 ES 



1 + 



(35) 



1 -r)P 

where P is the average of the day and the night survival probabilities. Since r ~ 0.16 and 
P « 0.3 for the best-fit solution, we see that the term in brackets in Eq. ( |35"D is about 
1.6. For the LMA solution, the best-fit predicted value for A is 8% for a 5 MeV recoil 
energy threshold, which corresponds to about 13% for SNO, in good agreement with the 
value of 12.5% given in Table |X[ There are small corrections to Eq.( [JSp due to the energy 
dependence of the various neutrino quantities (and the different locations on the earth of 
the SNO and the Super-Kamiokande detectors). 



C. The NC day-night effect 

There is no day-night effect in the NC for oscillations into active neutrinos. All active 
neutrinos are recorded with equal probability by the neutral current detectors. However, for 
oscillations into sterile neutrinos there can be a day-night effect since the daughter (sterile) 
neutrinos are not detectable. Thus a day-night effect in the NC would be a 'smoking gun' 
indication of sterile neutrino oscillations. 

For the region that is allowed at 99% CL by a global fit of the MSW sterile neutrino 
solution to all the available neutrino data Q, we find a NC MSW sterile neutrino day-night 
effect of 
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FIG. 9. The percentage difference between the night and the day CC rates. The figure shows 
for different neutrino scenarios the percentage difference between the predicted CC rate in SNO 
at night and the CC rate in the day [see Eq. (f28|)]. The error bars represent the 99% C.L. for 
the allowed regions of the six currently favored neutrino oscillation solutions Q. The results were 
calculated with a CC threshold of 5 MeV for the total electron recoil energy. 



A(NC, MSW Sterile) 



n nm +0.006 

u.uui_ 002 . 



(36) 



Although the predicted effect is small, it is important in principle since there are very few 
ways that sterile neutrino oscillations can be identified uniquely |26f. 



The neutral current day-night effect for solar neutrinos was first pointed out in Ref. [44 



Here we have calculated accurately the predicted range of the day-night asymmetry given 
the latest solar neutrino data, solar model, and a realistic model of the earth. 



X. SEASONAL EFFECTS 

We discuss in this section the seasonal dependences that are predicted by the currently 
favored neutrino oscillation solutions. We define a Winter-Summer Asymmetry by analogue 
with the Night-Day difference. Thus 

[Winter — Summer] 

%-s = 2 .— — -. (37) 

[Winter + Summer] 
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The earth's motion around the sun causes a seasonal dependence that can be calculated 
and is 



^4w-s, orbital = 0.064 (45 day averages) (38) 

for a 45 day Winter interval centered around December 21 and a 45 day Summer interval 
centered around June 21. The average length of the winter (summer) night during this 45 
day period is 15.4 (8.5) hours. The amplitude is reduced if the entire year is divided into 
two parts, with the winter average being taken as 182 days centered on December 21 and 
with the average length of the winter (summer) night being 14.1 (9.7) hours. In this case, 
the asymmetry is reduced by a factor of 1.5 from the 45-day value. Thus 

Av-s, orbital = 0.042 (182 day averages). (39) 

In what follows, we have removed the seasonal dependence due to the orbital motion 
from the quoted values of the seasonal dependence due to neutrino oscillation effects. 

Figure |TU| shows the predicted dependence upon the day of the year of the CC event rate, 
R cc , in SNO for each of the currently favored best-fit oscillation solutions. The vertical scales 
are different, reflecting the fact that the predicted seasonal variations are, e.g., relatively 
large for the best-fit LMA and VACl solutions, but are tiny for the MSW sterile solution. 



The annual average of the events rates shown in Fig. 10 yields the numbers shown in the 



second column of Table |vl[ The alert reader may notice that the phases of the variations in 
the two panels referring to vacuum oscillations are shifted by about two weeks with respect 
to the four panels that refer to MSW oscillations. This shift results from the fact that the 
earth and the sun are closest (relevant for vacuum oscillations) on January 4 and the day 
with the longest night is December 21 (relevant for MSW oscillations). 



Table IXII and Fig. ITTJ show the calculated percentage amplitudes for the 45 day winter- 



summer difference due to oscillations, A\y-s- In all cases, the best-fit oscillation solutions 
predict a winter-summer difference due to neutrino properties that is less than the orbital 
effects given in Eq. (|3~8|) and Eq. (|39|). Only rather extreme cases give amplitudes of Aw-s 
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FIG. 10. The seasonal dependence of the CC event rate in SNO. The figure shows the dependence 

on the time of the year of the predicted CC event rate in SNO for the six best-fit neutrino oscillation 

solutions described in Table |. The time labeled zero on the horizontal axis corresponds to January 

1. 

due to oscillations that are as large as the orbital amplitude, which will itself require a 
number of years to establish definitively p|. 



Table |XI] also gives the predicted values of ^4\v-s for a longer average, 182 of winter and 
182 of summer. For this case the statistical error will be reduced by about a factor of two, 
but the size of the effect is typically reduced by a factor of order 1.5 to 1.7. 

For the LMA solution, we showed in Ref. jKJ that to a good approximation A w _ s and 
A N _ D are related by the equation 



A 



w-s 



.4 



N-D 



24 h 



(40) 



where tw arid ts are the average lengths of the nights during the selected winter and summer 
periods, respectively. For the 45 day (182 day) intervals we are discussing here, the length 
at SNO of the winter night is 15.43 hours (14.10 hours) and the length of the summer night 
is 8.49 hours (9.72 hours). The term in brackets in Eq. (f40|) is 1.6 times larger for the 45 
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TABLE XI. The Winter-Summer Predicted Differ- 
ence in the CC Rate. The table gives the percentage 
difference, ^4\v-S! i n the Winter-Summer CC rates, defined 
by Eq. ([37]). The results on the first row for each oscillation 
solutions have been computed for 45 days in winter and 45 
days in the summer; the results on the second row are aver- 
age over 182 days of winter and 182 days of summer. The 
second column gives the best-fit value, and the third and four 
columns give the minimum and maximum values, allowed at 
99% CL for different neutrino oscillation solutions (cf. Fig. 
1 of Ref. Q). The winter-summer asymmetry due to the 
earth's motion around the sun, Ayv-S, orbital = 6.4%(4.2%) 
for a 45 day (182 day) average, has been removed from the 
values given here. 



Scenario 


b.f. 


min 


max 


Interval 


45 d 


45 d 


45 d 




182 d 


182 d 


182 d 


LMA 


3.6 


0.2 


7.7 




2.4 


0.1 


4.7 


SMA 


1.9 


-0.9 


11.8 




1.1 


-0.7 


7.8 


LOW 


1.0 


0.25 


2.9 




0.6 


0.15 


1.7 


VAC S 


3.4 


-0.9 


6.9 




2.2 


-0.6 


4.5 


VAC L 


2.0 


-0.6 


3.6 




1.3 


-0.4 


2.4 


MSW, Sterile 


-0.1 


-0.55 


1.1 




-0.1 


-0.3 


0.55 



day period (longer nights) than for the 182 day period. This accounts well for the ratios 
of Aw-s for the 45 day and the 182 periods that are given in Table [X| . Eq. (|40|) also 
produces well the individual values of Ayv-s for the LMA solution. Using the best-fit value 
of A N _ D = 12.48% (from Table [Xp and {t w - t s )/{24 h) = 0.29, we estimate A W -s = 3.6% 
for the 45 day average, in good agreement with the result given in Table pHf, 



tFor the LOW solution, Eq. (gOj) also gives a crude estimate of A\y-s, accurate to ~ 40%. The 
value of 8m? is smaller for the LOW than for the LMA solutions and therefore the typical oscillation 
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FIG. 11. The difference between the winter and the summer rates. The figure shows for different 
neutrino scenarios the percentage difference between the predicted CC rate in SNO for a 45 day 
interval during winter and the CC rate for a 45 day interval in the summer [see Eq. (37)]. The 
point labeled 'Orbital' represents the 45 day winter-summer difference due to the earth's motion 
around the sun; the amplitude of this orbital motion has been removed from the neutrino oscillation 
points. The error bars represent the 99% C.L. for the allowed regions of the six currently favored 
neutrino oscillation solutions The results were calculated with a CC threshold of 5 MeV for 
the total electron recoil energy. 



From the size of the predicted effects shown in Table X] and Fig. |TT|, we conclude that 
it will require many years of SNO operation to measure an accurate value of iw-s if the 
currently favored oscillation solutions are correct. 



XI. SMOKING GUN VS. SMOKING GUN 

What do we gain by combining the measurements of different smoking gun quantities? 
Once SNO has begun to report results for a variety of different quantities and an accurate 



length in matter is larger. The averaging of the oscillation effects required for the validity of Eq. (40) 



(see Ref. 15 ) is not complete for the LOW solution. 
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Monte Carlo of the experimental facility exists, then it will be possible to analyze simulta- 
neously a variety of different measurements using a global analysis method like Maximum 
Likelihood. In the meantime, we begin an initial illustrative exploration by analyzing pairs 
of SNO measurements. 

We show in this section how comparisons of the measurements of different smoking- 
gun quantities versus each other can enhance the deviation of a single measurement from 
the no-oscillation expectation and also shrink the globally-allowed range of the oscillation 
parameters. We concentrate on the most powerful pairwise combinations of variables. We do 
not illustrate all possible combinations, omitting some examples (like day-night effect versus 
first moment of the CC energy spectrum) that turn out to be less useful when examined 
quantitatively. 



We begin by displaying and discussing in Sec. PCI A| the predicted oscillation regions in 
planes defined by the double ratio [NC]/[CC] versus either 1) the day-night effect, A; 2) the 
first moment, (T), of the CC recoil energy spectrum; and 3) the neutrino-electron scattering 
reduced rate, [ES]. The double ratios involving the neutral current discriminate sharply 
between oscillation and no-oscillation scenarios and also reduce the range of acceptable 
oscillation parameters. In Sec. [XI B| and Sec. (XXCj, we discuss the location of the favored 
oscillation solutions in the [ES]/[CC] versus (T) plane and in the plane of the CC rate, Rcc 
versus the first moment, (T). 

For each plane defined by two SNO parameters and for each of the six neutrino oscillation 
solutions, we plot error bars that represent separately the 99% C.L. acceptable range of the 
neutrino parameters in the global fits to all the currently available solar neutrino data [|J. 
The la experimental uncertainties are summarized in Table |T|. The statistical uncertainties 
are computed assuming 5000 CC events, 1219 NC events, and 458 ES events. We assume 
that the hep uncertainties are symmetric and equal to the upper limit uncertainty, which 
slightly increases the error contours. For the no-oscillation case, only the experimental mea- 
surements are correlated. When neutrino oscillations occur, the predicted values for different 
parameters are also correlated. We include here only the correlations of the uncertainties 

45 



for the no-oscillation case; we do not include the correlated contours for the six different 
predicted oscillation solutions. A full calculation that includes the theoretical correlations 
between the different measured parameters and also includes asymmetric hep uncertainties 
should be carried out in the future, but this study is beyond the scope of the present paper. 

The correlations between different estimated experimental uncertainties cause the no- 
oscillation error ellipses to be tilted in Fig. |l2|-Fig. |16| . For purposes of illustration, we 
have assumed that the error correlations are as estimated in Ref. |Ll| . As we shall see from 
Fig. |T!-Fig. U, the tilt of the error ellipses can significantly influence the total statistical 
CL. assigned to a given set of results and therefore accurate determinations of the error 
correlations for the SNO experiment will be important. 



A. [NC]/[CC] double ratio versus other smoking guns 

1. [NCJ/[CCJ versus the day-night effect 

Figure [12] shows the values predicted by the different oscillation solutions in the plane 
of the [NC]/[CC] double ratio and the day-night asymmetry, An-d- Specifically, we plot 
the fractional shift in percent of the [NC]/[CC] double ratio from the standard model value 
of 0% on the vertical axis and the predicted value in percent of the day-night asymmetry 
An-d (standard model value of 0%) on the horizontal plane. Each of the currently allowed 
solutions, with the exception of the MSW Sterile solution, predicts points in the 5[NC]/[CC]- 
An-d plane that are more than 5a separated from the standard model solution (which is 
located at 0%, 0%). Moreover, the vacuum solutions are separated from the MSW solutions 
by amounts that exceed the expected measuring errors in An-d and [NC]/[CC]. It will, 
however, be more difficult to distinguish between different MSW solutions in the 5[NC]/[CC]- 
An-d plane, although some measured values would point to a unique solution. For example, 
a large positive value of A N _ D (> 20%) combined with a large value of [NC]/[CC] (> 2.3) 
would favor the LMA solution. The allowed region for the MSW Sterile solution is all 
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contained within the ellipse corresponding to the estimated 3a experimental uncertainty. 




20 40 

A N _ D (%) 

FIG. 12. The [NC]/[CC] double ratio versus the day-night asymmetry. The figure shows the 
fractional difference in percent of the neutral current versus charged current double ratio, [NC]/[CC] 
(Eq. |23|), from the no-oscillation value of 1.0, versus the day-night difference, ^n-d (Eq. p8|)in 
percent. The standard model value of ^4n-d is 0.0. For the six currently preferred oscillation 
solutions, the error bars represent the 99% C.L. regions for acceptable fits to all the available 
neutrino data. Estimated error contours for SNO are shown at the la, 2a, 3a, and 5a levels relative 
to the no-oscillation solution which lies at (0,0). 



2. [NC]/[CC] versus (T) 

Figure [13] shows the predictions of the different oscillation solutions in the 5[NC]/[CC] 
versus ST plane. All of the currently favored oscillation solutions, with the exception of the 
MSW sterile solution, predict locations in the 5[NC]/[CC] versus 5T plane that are sepa- 
rated by more than 5a from the standard model solution, which lies at (0.0, 0.0). However, 
the discrimination is almost entirely due to the [NC]/[CC] double ratio. The value of 5T 
only adds a large discrimination for the extreme VACs solution. The predicted values for 
the MSW Sterile solution extend out to 2.9%, which because of the correlation of the exper- 



47 



imental errors (which gives rise to the tilt of the error ellipses in Fig. [13]), can correspond to 
deviations as large as 3a from the no-oscillation solution. 




-5 5 

<5T (%) 

FIG. 13. The [NC]/[CC] double ratio versus the average recoil energy. The figure shows the 
fractional difference in percent of the neutral current versus charged current double ratio, [NC]/[CC] 
[Eq. (p3D], from the no-oscillation value of 1.0, versus the fractional difference in percent of the 
average electron recoil energy, (T) [Eq. (|TJ)] , from the no-oscillation value [Eq. flli])]. Contours 
at la, 2a, 3a, and 5a are shown relative to the no-oscillation solution at (0,0). 



3. [NC]/[CC] versus [ES] 

The most likely value for SNO to observe for the neutrino-electron scattering ratio [ES] 
is close to the Super-Kamiokande || value of [ES 5ft: ] = 0.475 (see for example Table |V11| or 
Fig. H). It is therefore convenient to define the quantity 5([ES]) as follows: 

[ES]o BS - 0.475 
v ; 0.475 V ; 

We have used parentheses rather than squared brackets in defining 5 (ES) because the shift 
in [ES] is measured relative to 0.475 rather than 1.0 . 



48 



!t 



400 



u 

£J 200 



LMA 



SMA 



LOW 




VAC 



VAC, 




MSW Sterile 

i i i L 



-20 20 

<5(ES) (%) 

FIG. 14. The [NC]/[CC] double ratio versus the neutrino-electron scattering double ratio. The 
figure shows the fractional difference in percent of the neutral current versus charged current double 
ratio, [NC]/[CC], from the no oscillation value of 1.0, versus the fractional shift in percent of the 
reduced neutrino-electron scattering ratio from the Super-Kamiokande value of 0.475 [see Eq. flUl)]. 
Contours at la, 2a, 3a, and 5a are shown relative to the no-oscillation solution at (0,0). 

Figure [14] shows the predictions of the different oscillation solutions in the 5[NC]/[CC] 



versus 5(ES) plane. Just as for Fig. [12| and Fig. [13|, all of the currently favored oscillation 
solutions, with the exception of the MSW Sterile solution, are well separated (more than 
5a away) from the standard model solution, which lies at (0.0,0.0). For some of the VACs 
solutions, the predicted large positive value of 5[ES] is incompatible with, and hence distin- 
guishable from, the predictions of the other currently allowed solutions. This discrimination 
is a result of combining the values of both [NC]/[CC] and [ES] since the measurement of 
either of these parameters by itself would not permit, according to Fig. the isolation of 
these VACs solutions. 
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FIG. 15. The [ES]/[CC] double ratio versus the average CC electron recoil energy. The figure 
shows the fractional difference in percent of the electron-scattering-CC double ratio, [ES]/[CC], 
from the no-oscillation solution of [ES]/[CC] = 1.0, versus the fractional shift in percent of the 
average electron recoil energy, (T) (Sec. pl[ ), from the no-oscillation value of Eq. (|l~4|). Contours 
at la, 2a, 3a, and 5a are shown relative to the no-oscillation solution at (0,0). 



B. Electron-scattering and CC double ratio versus CC energy spectrum 



Figure [15] shows the predictions of the different oscillation solutions in the 5[ES]/[CC] 
versus 5T plane. Although there are some predictions that extend well beyond the 5a 



contour in Fig. [15], these outlying predictions occur mostly for large values of [ES]/[CC] 
and should show up directly by comparing the neutrino-electron scattering rate with the 
CC (neutrino absorption) rate (see the discussion in Ref. [|J). The additional measurement 
of the first moment of the CC distribution, (T), does not add much to the discriminatory 
power of [ES]/[CC]. 
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FIG. 16. The CC rate versus the average recoil energy, (T). The figure shows the percent 
difference of the CC double rate, R cc , discussed extensively in Ref. 0, from the no-oscillation 
solution of R cc = 0.475, versus the fractional difference [see Eq. (|l|)] of the average electron recoil 
energy, (T), from the no-oscillation value of Eq. QT4|). Contours at la, 2a, 3a, and 5a are shown 
relative to the no-oscillation solution at (0,0). 

C. CC rate versus CC energy spectrum 

Figure [H] displays the six currently favored oscillation solutions in the plane of the CC 
rate, -Rcc ; and the first moment of the CC electron-recoil energy spectrum, (T). Some of 
the currently allowed VACs, SMA, and LMA solutions lie in this plane more than 5<r from 
the no-oscillation position. However, there are also currently allowed oscillation solutions 
that fall considerably closer to the no-oscillation point at (0.0,0.0). 



XII. SUMMARY AND DISCUSSION 

We concentrate in this section on describing the results for the predictions of the six 
currently favored neutrino oscillation solutions that are globally consistent at the 99% CL. 
with all of the available neutrino data. The neutrino solutions are described in Table [I] and 
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Fig. |. 

We begin this section by summarizing the results for parameters for which the estimated 
uncertainties are relatively small: 1) the neutral-current over charged current double ratio, 
Sec. pCll A| ; 2) the shape of the CC electron recoil energy spectrum, Sec. pCll B| ; 3) the day- 
night difference for the CC and for the NC, Sec. pCll C ; and 4) seasonal effects, Sec. |Xil D| . 
Altogether, we discuss six measurable quantities in Sec. pCII A| -Sec. pCII D| . 

We summarize the principal uncertainties, theoretical and experimental, in Sec. |XIIE| . 
The uncertainties due to the hep flux and the neutrino interaction cross sections are empha- 
sized in this section; the estimates of the experimental uncertainties are very preliminary. 

We then describe the predicted values and the potential inferences from SNO measure- 
ments for the CC rate, Sec. PCH F| , for the NC rate, Sec. PCH C| , and for the neutrino-electron 



scattering rate, Sec. PCIIH] . Next we discuss in Sec. PCH I| the neutrino-electron scattering to 
CC double ratio, which has some of the same advantages as the neutral-current to charged 
current double ratio, but suffers from a relatively large uncertainty in the CC interaction 
cross section. Finally we summarize in Sec. pCll J| our initial exploration of combining the 
analysis of different smoking gun indicators of neutrino oscillations. 



A. Neutral current over charged current double ratio: [NC]/[CC] 

All five of the currently favored oscillation solutions with active neutrinos predict that 
the double ratio, [NC]/[CC], will be separated from the no-oscillation value of 1.0 by more 
than 9a, estimated non-statistical errors. The uncertainties due to the cross sections and 
to the solar model almost cancel out of the double ratio. The minimum predicted value for 
[NC]/[CC] is 1.24 and the maximum predicted value is 5.15, all for a 5 MeV CC threshold. 
The estimated la total non-statistical error is only ±0.026; the statistical error will be the 
largest uncertainty unless more than 5.5 x 10 3 NC events are detected. The sterile neutrino 
solution lies in a disjoint region of [NC]/[CC] from 0.92 to 0.99. 

The results are summarized in Fig. [I] and are given in more detail, for two different 
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thresholds of the CC electron recoil energy, in Table [VIIl] of Sec. |VII . 

The double ratio [NC]/[CC] is an ideal smoking gun indicator of oscillations into active 
neutrinos. 



B. The shape of the CC electron recoil energy spectrum: (T) and <x(T) 

The shape of the CC electron recoil energy spectrum can be characterized by the first 
moment, (T), and the standard deviation, c(T), of the electron kinetic energy, T. With 
precision measurements of the spectrum, special features of the recoil energy spectra may 
also be detectable (cf. Fig. f|). 

If there are no oscillations, the first moment is (T)o = 7.422 x (1 ± 0.013) MeV, where 
the estimated uncertainty includes both the measurement and the calculational uncertain- 
ties. The best-estimate predictions for the different oscillation solutions correspond to a 
fractional shift between 0.1% (LMA and LOW solutions) to 3.7% (VAC S solution). The 
largest predicted value of the shift is 7.5% (VACs solution). The shift in the first moment 
may be measurable for the SMA, VACs, VACl, and MSW sterile neutrinos, but will be too 
small for a definitive measurement if the LMA or LOW solutions are correct. On the other 
hand, only the VACs solutions predict that the measured deviation of (T) from (T) will 
exceed three standard deviations for as much as half of the currently allowed solution space. 

Figure [5] and Table |IV| (of Sec. |T|) show the predicted range of shifts in the first moment, 
(T), of the recoil energy spectrum. A measurement of (T) to a la accuracy of ~ 100 keV 
will significantly reduce the allowed solution space for neutrino oscillations, but may not 
uniquely favor one particular solution. 

The calculated no-oscillation value of the standard deviation of the CC recoil energy 
spectrum is a = (<x 2 )o /2 = 1.852(1 ±0.049) MeV. It will be difficult to measure the predicted 
shifts from o"o since the total spread in shifts given in Table |V| is from —29 keV to 199 keV, 
while the estimated calculational and non-statistical measurement uncertainties are ±86 
keV. 
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C. Day-night difference: An— d and An_d(NC) 



In the absence of neutrino oscillations, there is a purely geometrical day-night difference 
that we have defined as the No Oscillation (NO) effect and whose value for the SNO detector 
we have given in Eq. fl29|). We have removed the NO effect from all of the calculated day- 
night values given in this paper. 

For the currently favored MSW active neutrino solutions, the best-fit predictions for 
the average difference between the nighttime and the daytime CC rates, An-d, vary from 
2% for the SMA solution to 12.5% for the LMA solution, all for a 5 MeV recoil electron 
energy threshold. Small values (< 1%) of the day-night difference would be consistent with 
any of the three MSW solutions, but very large values of the day-night difference are only 
expected for some of the LMA solutions. The maximum expected difference for the most 
extreme LMA solution is 28.5%, whereas the maximum difference expected for the LOW 
(SMA) solution is 13.5% (12%). The MSW sterile solution predicts values for the day-night 
asymmetry between —0.5% and +1.1%. Table [X| presents similar results also for an 8 MeV 
electron recoil energy threshold. 

For vacuum oscillations, the predicted values of the day-night effect are small, but non- 
zero (see Table 0). 

Initially, the dominant uncertainty for the day-night effect will be purely statistical. The 
most difficult problems will ultimately arise from systematic effects, such as the symmetry 
of the detector and the separation of the CC events from NC and scattering events, that 
will have to be modeled by a detailed SNO Monte Carlo simulation. The purely statistical 
error may be of order 4% after one full year of operation. Whether or not the systematic 
errors affect in an important way the measurement of the day-night effect will depend upon 
the actual magnitude of An-d and the size of the systematic uncertainties. 

Figure ^ and Table |X| present the numerical results for the CC day-night asymmetry 



which is defined by Eq. ( pq) of Sec. |X[ If one of the MSW active neutrino solutions is 
correct, then the day-night difference could become apparent early in the operation of SNO. 
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This possibility exists for the MSW active solutions, but is not required. 

There is no day-night effect in the NC for oscillations into active neutrinos. Oscillations 
into sterile neutrinos give a small effect: A N _ D (NC, MSW Sterile) = - COOlt^ [see 
Eq. fl36Dl- This effect is important in principle, since it is a clear distinction between active 
neutrinos and sterile neutrinos. However, the predicted size is too small to be measured 
with SNO. 

All of the currently favored neutrino oscillation solutions (MSW or vacuum oscillations 
into active neutrinos, as well as MSW oscillations into sterile neutrinos), predict that 

| A N _ D (NC) | < O.Of. (42) 

We have written Eq. fl4~2"D in its most general form. Of course, v4n-d(NC) is predicted to be 
identically zero for all neutrino oscillations (vacuum or MSW) into active neutrinos. 

The measurement by SNO of An-d(NC) will be an important test of neutrino oscillation 
models. If we obtain independent evidence that solar neutrino oscillations involve active 
neutrinos, or if one hypothesizes that sterile neutrinos play no role in solar neutrino os- 
cillations, then the measurement of An-d(NC) can be regarded as a test of the standard 
electroweak model. 



D. Seasonal effects: Aw— s 



Figure [TTjand Table |XI]give the amplitudes of the CC winter-summer differences, Aw_s, 
that are predicted by the favored neutrino oscillation solutions. The results can be compared 
with the amplitudes expected from the orbital motion of the earth, which are given in 
Eq. (|38| ) and Eq. (^) . In all cases, the current best-fit oscillation solutions predict a winter- 
summer amplitude that is less than the amplitude due to the earth's orbital motion. Only for 
a small fraction of the currently allowed oscillation parameters does the predicted amplitude 
due to oscillations exceed the amplitude due to the earth's orbital motion. 

We conclude that it will probably be difficult to measure Aw-s- However, we note that 
the prediction that Aw-s is small is a prediction that can and should be tested. 
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E. Uncertainties 



Table [TTT] presents a convenient summary of the estimated calculational and measurement 
uncertainties for different experimental quantities that will be determined by SNO. We 
present in Sec. || a full description of how we estimate these uncertainties. We do not include 
the effects of background events; there is no reliable way of estimating the background prior 
to actual measurements in the SNO detector. We also do not include misclassification 
uncertainties, e.g., ES events mistaken for CC events or NC events mistaken for CC events. 
These errors must be determined by the detailed SNO Monte Carlo simulations. 

The quantitative influence of the hep flux of neutrinos on the measurement accuracy 
of different quantities has been evaluated here for the first time. In addition, we include 
estimates of uncertainties due to the width of the resolution function for the recoil electron 
energies, the absolute energy scale, the 8 B neutrino energy spectrum, the interaction cross 
sections, and the number of events counted. 

Our present limited experimental knowledge of the hep flux causes an uncertainty of 
~ 2% in all three of the rates that will be measured by SNO, i.e., the CC rate, the NC rate, 
and the neutrino-electron scattering rate (see Table However, measurements of the CC 
spectrum by SNO can reduce the uncertainty in the hep flux |3J and therefore decrease the 
contribution of the hep to the error budgets of different SNO measurables. Most recently, 



an improved theoretical calculation [f|6(] of the low energy cross section factor for the hep 
reaction has increased the best-estimate of the flux from the standard solar model by about 
a factor of five relative to the value given in BP98 and used in the present paper. 

The neutrino cross section uncertainties for the CC and the NC rates are the largest 



entries in Table pj. The lack of knowledge of these interaction cross sections will limit the 
interpretation of the measured rates of both the CC and the NC interactions. 

Table || summarizes the recent cross section calculations for the CC, NC, and for the 
ratio, NC/CC. We have computed uncertainties in measurable quantities using the entries 
in Table [TI| and the algorithm given in Eq. (0). It is a matter of judgment as to how many 
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standard deviations should be assigned to the difference computed from Eq. (0). We believe 
that we are being reasonable and conservative in regarding this difference as la. But, we 
stress the need to greatly increase the limited number of entries in Table [TI| so that a more 
informed estimate of the uncertainties can be made. The need for additional calculations 
is particularly urgent for the double reduced ratio [NC]/[CC]. We have used just the two 
entries in Table |T| to estimate that the uncertainty in [NC]/[CC] is an order of magnitude 
less than the separate uncertainties in [NC] and [CC] . While plausible, it is essential to check 
that this cancellation of uncertainties in the double ratio is indeed a general characteristic 
of accurate calculations of the neutrino interaction cross sections. 



Towner j|7]] has computed radiative corrections for both the CC and the NC neutrino 
reactions on deuterium. The effect of the radiative corrections is generally small. Radiative 
corrections change the first moment by about 0.1% and the second moment by about 0.3%. 
Although not computed by Towner, the effect of radiative corrections on time- dependent 
quantities such as the day-night effect, the zenith angle distribution, and the seasonal effects 
is expected to be similarly small (< 0.3%). For the double ratio [NC]/[CC], the effect of the 
radiative corrections is larger, 0.5%, if the photons from the inner bremsstrahlung in the 
CC reactions are not detected. In the extreme case in which all of the inner bremsstrahlung 
photons are somehow detected by SNO, the [NC]/[CC] ratio would be increased by 4% |. 

The knowledge of the cross section uncertainties can be improved by further calculations, 
especially those based upon chiral symmetry. Calculations should be carried out for a variety 
of models and approximations and with the full range of allowed nuclear and particle physics 
parameters. An initial step in this direction has been taken by Butler and Chen ||31|| , who 
have calculated the NC reaction in effective field theory. A full exploration of the allowed 



^For neutrino-electron scattering, radiative corrections are also small and have been computed 
explicitly by Bahcall, Kamionkowski, and Sirlin [fj^]. The Super-Kamiokande collaboration includes 
these calculated radiative corrections in their analyses ||. 
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range of CC and NC cross sections for neutrinos incident on deuterium is an urgent and 
important task for the theoretical nuclear physics community. 

Further experimental work on neutrino interactions with deuterium would be extremely 
valuable. Butler and Chen [|3T] have pointed out that a measurement of the two-body 
matrix element could be made using the reaction e + 2 H— > e + n + p. This is an urgent and 
important task for the experimental nuclear physics community. More precise measurements 
of the anti-neutrino disintegration cross sections made with reactors would be valuable in 
choosing between and guiding theoretical calculations (for a state-of-the-art discussion of 



the experimental possibilities see Ref. p2|). It would also be useful to test the accuracy 
of the calculational procedures, albeit at higher neutrino energies, by performing neutrino 
absorption and disassociation experiments on deuterium with a stopped muon beam. 

F. The CC rate: [CC] 

The charged current rate will be one of the first results to be obtained with SNO. The 
reduced neutrino-absorption rate, [CC], defined by Eq. ([T7|), can be compared with the 
neutrino-electron scattering ratio measured by Super-Kamiokande ||, 0.475 ± 0.015. If the 
CC ratio is measured to be less than 0.475, then that would be evidence that u-e scattering 
includes contributions from muon or tau neutrinos and therefore neutrino oscillations are 
occurring. 



Table |V| presents the predicted CC reduced rates, [CC], for the six currently favored 
oscillation solutions. For example, the predicted ratio for oscillations into active neutrinos 
ranges from 0.29 (LMA, minimum value) to 0.46 (SMA, maximum value), if the recoil 
electron energy threshold is set at 5 MeV. There is about an equal chance, according to 
Table |VT] and Fig. 2 of Ref. that the measured value of [CC] will lie more than 3a 
below the no-oscillation value of 0.475. Four of the solutions, the LMA, SMA, LOW, and 
VACs solutions, all have best-fit global solutions that predict [CC] < 0.40 and each of these 
solutions has some region of neutrino parameter space that gives values as low as 0.35 or 
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below. 

The MSW sterile solution predicts values for [CC] that are very close to 0.475 for an 
electron recoil energy threshold of 5 MeV. For a threshold of 8 MeV, the MSW sterile solution 
predicts values for [CC] that even exceed 0.475 (see explanation of this interesting fact in 
Ref. [|). 

The discriminatory power of the CC rate measurement could be increased significantly 
if the uncertainty in the CC cross section could be decreased (see discussion in Sec. [XII E| 
above) . 



G. The NC rate: [NC] 

The reduced neutral current rate, [NC], should be equal to 1.0 if the standard solar 
model prediction of the 8 B flux is exactly correct and if there are no oscillations into sterile 
neutrinos. Oscillations into active neutrinos would preserve the neutral current rate and 
would not change the 1.0 predicted value of the reduced rate. 

The interaction cross section constitutes the largest uncertainty, ±6%, in determining 
experimentally the reduced neutral current rate. The uncertainties in calculating the solar 
flux (+18%, —16%, see Ref. ||) provide the biggest complication in interpreting the neutral 



current measurement directly in terms of neutrino physics [see Eq. (|20|) for the uncertainties 
in measuring and interpreting [NC]]. 

The sensitivity of [NC] to the true solar flux is a problem for particle physics, but an 
advantage for astrophysics. 

The measurement of the neutral current rate will provide crucial information about 
the true 8 B solar neutrino flux provided there are no oscillations into sterile neutrinos. 
Unfortunately, the MSW sterile neutrino prediction, [NC] Ste riie = 0.48 ± 0.01 [see Eq. ([21])], 
is within about 3a of the no-oscillation value of 1.0 when one includes the uncertainty in 
the solar model flux. Hopefully, the uncertainty in the predicted value of the standard 
solar model flux will be reduced somewhat by precise laboratory measurements of the 8 B 
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production cross section that are now in progress. 

A measurement of [NC] larger than or close to 1.0 would be evidence against sterile 
neutrino oscillations and would support the solar model estimate for the 8 B flux (provided 
the experimental value is not larger than 1.5). A measurement between 1.0 and 0.5 could be 
interpreted as indicating a solar flux somewhat lower than the best estimate or as providing 
evidence for sterile neutrinos. A measurement significantly below 0.5 would be a clear 
indication of oscillations into sterile neutrinos, but would conflict with the Kamiokande f|| 
and Super-Kamiokande || measurements of the v-e scattering rate. 



Bilenky and Giunti [48[ have pointed out that a comparison of the time-dependence of 
the CC and NC SNO rates on the time scale of the llyear solar cycle could test the spin- 
flavor precession scenario. According to this hypothesis, the NC rate would remain constant 
throughout the solar cycle while the CC rate would vary with phase in the cycle. 

H. The neutrino-electron scattering rate: [ES] 

Table |V11| and Fig. |6| show the predicted values of the reduced neutrino-electron scattering 
ratio, [ES], for the six currently favored oscillation solutions. Not surprisingly, the predicted 
ratios cluster close to the Super-Kamiokande value of 0.48 ||. The most extreme values range 
from 0.45 (minimum allowed for the MSW sterile solution) and 0.52 (maximum allowed for 
the VACs solution). These values are all for a recoil electron energy threshold of 5 MeV 



(see Table |VII| for results for an 8 MeV threshold) and a 99% C.L. for the allowed range of 
oscillation solutions. 

For the first five or ten years of operation of SNO, the dominant measurement uncertainty 
for [ES] will be statistical (cf. Table |VTlD . The observed rate of neutrino-electron scattering 
events in SNO is expected to be only ~ 10% of the CC rate. 

Although the SNO detector is different from either the Kamiokande or the Super- 
Kamiokande detectors, and the value of [ES] depends somewhat on threshold (see Table |V11| ) 
and on the instrumental parameters such as the width of the energy resolution function, the 
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bottom-line results for [ES] should be similar in all cases for these water Cherenkov detectors. 

The SNO measurement of v-e scattering will provide an important check of SNO versus 
Kamiokande and Super-Kamiokande and vice-versa. In addition, the value of [ES] as deter- 
mined in SNO can be used in connection with other SNO measurements to constrain the 
allowed neutrino parameter oscillation space. 



I. The v-e to CC double ratio: [ES]/[CC] 

The double ratio of [ES]/[CC] is, like the [CC]/[NCC] double ratio, largely insensitive to 



solar model uncertainties (see Table |HI|) . Moreover, some of the systematic measurement 
uncertainties are reduced because the same techniques are used to detect v-e scatterings 
([ES]) and neutrino absorption ([CC]). The principal difficulty in interpreting measurements 
of the double ratio [ES]/[CC] at the present time is the large uncertainty, 5.8%, in the CC 
reaction cross section. This uncertainty is almost six times larger than any other known 
contributor to the [ES]/[CC] error budget [see Eq. (pTj)]. 

Figure || and Table [IX] present the predicted range of values for [ES]/[CC] for the six 
currently favored neutrino oscillation solutions. For oscillations into active neutrinos, 1.03 < 
[ES]/[CC] < 1.65 for a 5 MeV recoil electron energy threshold. The corresponding limits 



are 1.05 and 1.67 for an 8 MeV energy threshold (see Table |Xj). The total non-statistical 



uncertainty is estimated to be 7% [see Eq. (|7|)] and the statistical uncertainty will be about 
5% after the accumulation of 5000 CC events. 

We conclude that Nature has adequate opportunity to choose an oscillation solution into 
active neutrinos in which the ratio [ES]/[CC] is many sigma from the no-oscillation value 
of 1.0. Nevertheless, the contrast between the no-oscillation value and the currently favored 
oscillation predictions is much less for [ES]/[CC] than it is for [NC]/[CC]. The greater power 
of [NC]/[CC] can be seen most clearly by comparing Fig. [7], which has a vertical scale that 
extends from 0.5 to 7.0, with Fig. 0, which has a vertical scale that extends only to 3.0. 



61 



J. Smoking gun versus smoking gun 

The full diagnostic power of SNO will be achieved by analyzing simultaneously all of 
the measurements, including upper limits. This full analysis requires detailed and mature 
Monte Carlo simulations based upon experimental calibrations. 

Figure ^-Figure ^ provide an illustrative introduction of what can ultimately be 
achieved by the simultaneous analysis of the full set of SNO measurables. The figures are 
two dimensional slices in the multi-dimensional SNO parameter space; we plot one smok- 
ing gun against another smoking gun. Contours ranging from la to 5cr are shown for the 
no-oscillation case and include estimates for the error correlations. The error bars for the 
different oscillation scenarios represent the range of values predicted for each smoking gun 
independently. 

XIII. WHAT ARE OUR MOST IMPORTANT CONCLUSIONS? 

The paper contains many specific results. Here is our personal list of our most important 
conclusions. 

(1) The neutral current to charged current double ratio. All currently favored 
active neutrino oscillation solutions predict a value for the double ratio, [NC]/[CC], of neutral 
current to charged current event rates that is, with our best estimates for the theoretical and 
experimental uncertainties, more than 9a away from the no-oscillation solution (neglecting 
statistical uncertainties). If statistical uncertainties are included for 5000 CC events (and 
1219 NC events), then the minimum discrepancy is reduced to 6a. 

(2) Day-Night differences in the CC rate. Large differences are predicted between 
the day and the night CC rates for some currently favored MS W solutions. For a 5 MeV 
electron recoil energy threshold, the best-fit differences vary between —0.1% (MSW Sterile) 
and 12.5% (LMA). The largest predicted value among all the currently allowed solutions 
is 28.5% (LMA), which could be detectable in the first year of operation of SNO. Similar 
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results are predicted for an 8 MeV recoil energy threshold. Vacuum oscillations have average 
day-night differences of order 1% or even less. Small values (< 1%) of the CC day-night rate 
difference would be consistent with any of the six currently favored two-flavor oscillation 
solutions. The day-night difference of the NC is predicted to be < 0.01% for all solutions 
(and is non-zero for the MSW Sterile solution). 

(3) Uncertainties. The uncertainties in the absolute values of the neutrino cross sections 
for the CC and for the NC current reactions are the largest known uncertainties. These 
uncertainties limit the interpretation of the separate CC and NC rates, but cancel out (to 
an accuracy of better than 1%) of the [NC]/[CC] ratio. 

(4) Specrum distortion. The first moment of the CC electron recoil energy spectrum 
describes well the predicted deviations for all except the VACl solution. For all the MSW 
solutions and for the VACs solution, the predicted spectrum distortion is smooth and mono- 
tonic in the region accessible to SNO (see Fig. f|) and hence can be characterized by a single 
parameter. The currently favored oscillation solutions predict a range of deviations of the 
first moment, most of which are less than the estimated 3er experimental uncertainty. The 
largest predicted deviations are for the VACs (best-fit predicted deviation 283 keV) and 
SMA (best-fit predicted deviation 218 keV) solutions. 

The VACl solution generically predicts a bump and a dip in the low and middle energy 
region of the SNO electron recoil energy spectrum (see Fig. |j). For the best-fit VACl 
solution, this modulation is about 30% and if observed would be strong evidence for the 
VACl scenario. The predicted modulation occurs in a region where the event rate is expected 
to be relatively high. 

(5) Characteristic size of effects. The current best-fit global neutrino oscillation so- 
lutions typically predict small effects, of order several percent or less, for all of the quantities 
that are sensitive to oscillations which SNO will measure, except [NC]/[CC]. However, for 
some allowed oscillation solutions, the difference between the day and the night rates and 
the distortion of the shape of the CC electron recoil energy spectrum may be relatively large. 

(6) Sterile neutrinos. The current best-fit MSW Sterile solution predicts, relative 
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to the no-oscillation solution, a 1.8% shift in the first moment of the CC electron recoil 
energy spectrum and a CC rate that is larger than for the other currently allowed oscillation 
solutions. The neutral current rate is predicted to be 0.465 ± 0.01 of the standard solar 
model rate, i.e., the MSW Sterile solution predicts a much smaller value for the neutral 
current rate than the other allowed oscillation solutions. The sterile solution also predicts 
a small but non-zero value for the difference between the NC rate during the day and the 
NC rate at night. The CC day-night difference is predicted to be small, but not as small as 
for the NC day-night difference. It will be difficult to discriminate with SNO between the 
no-oscillation solution and the currently allowed MSW Sterile solution. 

The Sudbury Neutrino Observatory will enrich particle physics with measurements of 
many effects that are sensitive, in different ways, to neutrino oscillations . 
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addresses. We are indebted to E. Akhmedov, E. Beier, S. Bilenky, D. Cowan, E. Kearns, J. 
Feng, M. Fukugita, K. Kubodera, E. Lisi, A. McDonald, and Y. Nir for valuable comments 
on a draft copy of this manuscript. JNB and AYS acknowledge partial support from NSF 
grant No. PHY95-13835 to the Institute for Advanced Study and PIK acknowledges support 
from NSF grant No. PHY95-13835 and NSF grant No. PHY-9605140. 
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